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We upgrade a SUe quark-model description for the nucleon-nucleon and hyperon-nucleon inter- 
actions by improving the effective meson-exchange potentials acting between quarks. For the scalar- 
and vector-meson exchanges, the momentum-dependent higher-order term is incorporated to reduce 
the attractive effect of the central interaction at higher energies. The single-particle potentials of 
the nucleon and A, predicted by the G-matrix calculation, now have proper repulsive behavior in 
the momentum region gi = 5 - 20 fm~^. A moderate contribution of the spin-orbit interaction from 
the scalar-meson exchange is also included. As to the vector mesons, a dominant contribution is the 
quadratic spin-orbit force generated from the p-meson exchange. The nucleon-nucleon phase shifts 
^ ■ at the non-relativistic energies up to Tiab = 350 MeV are greatly improved especially for the 

states. The low-energy observables of the nucleon-nucleon and the hyperon-nucleon interactions are 
also reexamined. The isospin symmetry breaking and the Coulomb effect are properly incorporated 
in the particle basis. The essential feature of the AN-T,N coupling is qualitatively similar to that 
obtained from the previous models. The nuclear saturation properties and the single-particle poten- 
tials of the nucleon. A, and E are reexamined through the G-matrix calculation. The single-particle 
' potential of the E hyperon is weakly repulsive in symmetric nuclear matter. The single-particle 

spin-orbit strength for the A particle is very small, in comparison with that of the nucleons, due to 
■ the strong antisymmetric spin-orbit force generated from the Fermi-Breit interaction. 
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I. INTRODUCTION 

Q 

^3 , . 

One of the most important purposes of studying the nucleon-nucleon (NN) and hyperon-nucleon (YN) interactions 
• • in the quark model is to understand comprehensively the fundamental strong interaction in a natural picture, in which 
the quark-gluon degree of freedom is relevant to describe the short-range part of the interaction, while the medium- 
and long-range parts of the interaction are dominated by the meson-exchange processes. We have recently achieved a 
■ simultaneous and realistic description of the NN and YN interactions in the resonating-group (RGM) formalism of the 
spin-flavor SUq quark model In this approach the effective quark-quark (qq) interaction is built by combining a 

phenomenological quark-confining potential and the colored version of the Fermi-Breit (FB) interaction with minimum 
effective meson-exchange potentials (EMEP) of scalar and pseudo-scalar meson nonets directly coupled to quarks. 
The flavor symmetry breaking for the YN system is explicitly introduced through the quark-mass dependence of the 
Hamiltonian. An advantage of introducing the EMEP at the quark level lies in the stringent relationship of the flavor 
dependence appearing in the various NN and YN interaction pieces. In this way we can utilize our rich knowledge of 
the NN interaction to minimize the ambiguity of model parameters, which is crucial since the present experimental 
data for the YN interaction are still very scarce. 

In this study we upgrade our model ||l|-g by incorporating such interaction pieces provided by scalar and vector 
mesons as the spin-orbit (LS), quadratic spin-orbit {QLS), and the momentum-dependent Bryan-Scott terms. Intro- 
duction of these pieces to the EMEP is primarily motivated by the insufficient description of the experimental data 
by previous models. First, some discrepancy of the NN phase shifts in previous models requires the introduction 
of vector mesons. For example, phase shift in the model ESS 0] is more attractive than experiment by 10° 
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around Tiab ''^ 300 MeV. This implies that the one-pion tensor force is too strong in our previous models. In the 
standard one-boson exchange potentials (OBEP), the strong one-pion tensor force is partially weakened by the p 
meson tensor force. We use the QLS force of vector mesons from the reasons given below. Furthermore, some phase 
shifts of other partial waves deviate from the empirical ones by a couple of degrees. Another improvement is required 
as for the central attraction. The G-matrix calculation using the quark-exchange kernel explicitly |^ shows that 
energy-independent attraction, dominated by e-meson exchange, is unrealistic, since in our previous models the single 
particle (s.p.) potentials in symmetric nuclear matter show a strongly attractive behavior in the momentum region 
gi = 5 - 20 fm~^. We have shown in |Q that this flaw can be removed by introducing the momentum-dependent 
higher-order term of scalar-meson exchange potentials, the importance of which was first pointed out by Bryan and 
Scott 1^. In the higher energy region, the LS term of the scalar mesons also makes an appreciable contribution, in 
addition to this momentum dependent term. 

Another purpose of the present investigation is to examine the charge symmetry breaking (CSB) and the Coulomb 
effect from the viewpoint of the quark model. It is well known that the ^5*0 phase shift of the pp interaction is slightly 
less attractive than that of the np interaction. This charge independence breaking (CIB) is partially explained by the 
so-called pion-Coulomb correction [0, which implies 1) the small mass difference of the neutron and the proton, 2) 
the mass difference of the charged pion and the neutral pion, and 3) the Coulomb effect. Furthermore, it was claimed 
long ago that the Ap interaction should be more attractive than the An interaction, since the binding energy of the 0+ 
ground state of ^Hc is fairly larger than that of The CSB energy of 350 keV in these isodoublet hypernuclei 

is much larger than ~ 100 keV CSB effect seen in the '^H-'^IIe binding energy difference after the correction of the 
pp Coulomb energy in ■^He is made. The early version of the Nijmegen potential already focused on this CSB in 
the OBEP including the pion-Coulomb correction and the correct threshold energies of the AN-'EN coupling in the 
particle basis. The RGM calculation using the particle basis is rather cumbersome, since all the spin-flavor factors of 
the quark-exchange kernel should be recalculated by properly incorporating the z-components of the isospin quantum 
numbers. Furthermore, there is a problem inherent in the RGM formalism: the internal energies of the clusters 
are usually not properly reproduced when a unique model Hamiltonian is used. We have given in jl^ a convenient 
prescription to avoid this problem without spoiling the exact effect of the Pauli principle. For the Coulomb effect, 
we calculate the full exchange kernel without any approximation. The pion-Coulomb correction and the correct 
treatment of the threshold energies in the particle basis are found to be very important for the detailed description 
of the low-energy observables in the T,N - AN coupled-channel problem. 

With these renovations of EMEP and the framework, we have again searched for model parameters in the isospin 
basis to fit the most recent result of the NN phase shifts, the deuteron binding energy, the ^5*0 NN scattering length, 
and the low-energy YN total cross section data. This model is named fss2 since it is based on our previous model 
FSS |^-||]. The agreement of the phase-shift parameters in the NN sector is greatly improved. The model fss2 shares 
the good reproduction of the YN scattering data and the essential features of the AN-Y,N coupling with our previous 
models The single-particle (s.p.) potentials of N, A and S are predicted through the G-matrix calculation. |^ 

The strength of the s.p. spin-orbit potential is also examined by using these G-matrices [ |lO| |. 

In the next section we first recapitulate the formulation of the {3q)-{3q) Lippmann-Schwinger RGM (LS-RGM) 
and the G-matrix calculation Q using the explicit quark-exchange kernel. Section II B introduces a new EMEP 
Hamiltonian for fss2 in the momentum representation. This serves to clarify the difference between the present model 
fss2 and the previous two models, FSS and RGM-H ||^||]. The spatial part of the quark-exchange kernel and the spin- 
flavor factors in the EMEP sector are given in Appendices A and B, respectively. The model parameters determined in 
the isospin basis are discussed in Sec. II C. Short comments are given in Sec. II D with respect to the special treatment 
in the particle basis, including the Coulomb force in the momentum representation. Section III presents results and 
discussions. We first discuss in Sec. Ill A the NN phase shifts, differential cross sections and the polarization for the 
energies Tiab < 800 MeV. Special attention is paid to the effect of inelastic channels, which is not taken into account 
in the present framework. The five invariant amplitudes for the pp scattering are also examined at the highest energy 
Tiab = 800 MeV, in order to clarify the behavior of the s.p. potentials in the asymptotic momentum region and to find 
a clue to the missing ingredients in the present framework. The deuteron properties and the effective-range parameters 
of the NN system are discussed in Sec. Ill B. A simple parameterization of the deuteron wave functions is given in 
Appendix C. Sections III C, D, and E discuss the phase-shift behavior and the characteristic features of the S^p, AN, 
and S~p systems, respectively. In Sec. HIE, the pion-Coulomb correction and the inelastic capture ratios at 
rest and in flight are also discussed in the particle basis. The YN cross sections in the low- and intermediate-energy 
regions are discussed in Sec. Ill F. The G-matrix calculation using fss2 is presented in Sec. Ill G. This includes the 
discussion of the nuclear saturation curve, the density dependence of the s.p. potentials and the Scheerbaum factors 
of the s.p. spin-orbit strength in symmetric nuclear matter. The flnal section is devoted to a summary. 
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II. FORMULATION 



A. The Lippmann-Schwinger formalism for (3q) - (3q) RGM and the G-matrix equation 




A new version of our quark model employs the Hamiltonian which includes the interactions generated from the 
scalar (S), pseudoscalar (PS) and vector (V) meson exchange potentials acting between quarks: 

H = j:{rn.c^ + ^-Ta]+X.\U^+U:;r + l^U^ + l^U^^' + l^U:^'\- (2-1) 

Here Up/ is a confinement potential with a quadratic power law, and Ul/ is the full FB interaction with explicit 
quark-mass dependence. It is important to note that this confinement potential gives a vanishing contribution to the 
baryon-baryon interaction, since we assume {Os)^ harmonic oscillator wave functions with a common width parameter 
b for the internal cluster wave functions. Also, all the contributions from the FB interaction are generated from the 
quark-exchange diagrams, since we assume color-singlet cluster wave functions. These features are all explained in our 
previous publications Q|. When the calculations are made in the particle basis, the Coulomb force is also introduced 
at the quark level. The RGM equation for the parity-projected relative wave function XaiR) derived from the 
variational principle {S'^\E ~ H\'i>) — 0, and it reads 



(liR) - ^IrfJ?' Gac'iRR^E) xl'iR) , (2.2) 



where Gaa' {R, R ', E) is composed of various pieces of the interaction kernels as well as the direct potentials of EMEP: 

G„„, {R, R-E)^ 5{R ~R)J2T. ^SL' (R) + E - ■ (2.3) 

13 n o 

The subscript a stands for a set of quantum numbers of the channel wave function; a ~ [1/2(11) ai, l/2(ll)a2] 
SSzYIIz;V, where l/2(ll)a is the spin and SU3 quantum number in the Elliott notation (A/x), a {— YI) is the flavor 
label of the octet baryons {N = 1(1/2), A — 00, E = 01 and S — —1(1/2)), and V is the flavor-exchange phase 
In the NN system with 0102 = NN, V becomes redundant, since it is uniquely determined by the isospin 
as V = (— 1)^"'^. These are the channel specification scheme in the isospin basis. In the particle basis, necessary 
modification should be made for the flavor degree of freedom. The relative energy in the channel a is related to 
the total energy E of the system in the center-of-mass (cm.) system through Sa = E — E™^. Here E™*^ — E™*^ + E™^ 



with a = aia2- In Eq. (2_^) the sum over Q for the direct term implies various contributions of interaction types for 
the meson-exchange potentials, while (3 specifies the meson species. On the other hand, 17 for the exchange kernel 
A42^,{R, R) involves not only the exchange kinetic-energy ( K) te rm but also various pieces of the FB interaction, 
as well as several components of EMEP. The RGM equation ( |2.2| ) is solved in the Lippmann-Schwinger formalism 
developed in Q (which we call LS-RGM). In this formalism, we first calculate the basic Born kernel defined through 

M^AQf,Q^■,E) = {e''>r^\Go.„iR,R;E)\e^'^-''') 

= E E ^Caa'(g/, 90 + E 90 O'\qf,q0 - M^^,{qf,q,) , (2.4) 



where is the relative energy in the final channel (in the prior form). Each component of the Born kernel Eq. (2.4) 
is given in terms of the transferred momentum k — qf ~ qt and the local momentum q = {qj + qi)/2. In Eq. (2.4) 
the space-spin invariants — 0^{qf, qi) are given by = 1 and 

with n=[q,xqf], S= i(^i + ^2) , S^--^ = i(^i - ^2) , Pa = ^ ^ ' • (2.5) 
For the tensor and QLS parts, it would be convenient to take four natural operators defined by 

= 5i2(fc, fe) , O^' = Si2{q, q) , O^" = S^ik, q) , O'?^'^ = 5i2(n, n) , (2.6) 
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where Si2{a,b) = (3/2)[(cri • a)((T2 • b) + {ct2 ■ a){(Ti • b) ] - (cri • <T2)(a • b). The direct Born kernel M^1^^,{qf , q.^) in 
Eq. (2^) is exphcitly given in Sec. IIB. The exchange Born kernel Mf^,{qf,qi) is given in Appendix B of for the 



FB interaction and in Appendix A for the EMEP. The LS-RGM equation is given by 

kl-k"^ + is 



T^a{p,q;E) = V^aip,q;E)+J2j^^ J dk V^fs{p,k; E)- 



Tf3o.{k,q;E) 



where the quasi-potential Vya{p, Q', E) or more generally V^p{p, q; E) is calculated from 
V,p{k,k'-E) = \ [ M^p{k,k'-E) + {-iff^Vp M^f,{k,-k'-E) ] . 



(2.7) 



(2.8) 



After the standard procedure of the partial- wave decomposition]^ the LS-RGM equation (2/7) is solved by the Noyes- 
Kowalski method [^,0. The singularity at fc = fc^ is avoided by separating the momentum region into two pieces. 
The intermediate fc-integral over < fc < fc^ is carried out using the Gauss-Legendre 15-point quadrature formula 
and the integral over kp < k < oo using the Gauss-Legendre 30-point quadrature formula through the mapping 
k = kp + tan(7r(l + x)/4). 



The LS-RGM equation (2.7) is straightforwardly extended to the G-matrix equation by a trivial replacement of the 



free propergator with the ratio of the angle-averaged Pauli operator and the energy denominator: 
G^^[p,q;K,u;) = V^^{p,q-E) + Y.Tfyz I ^ k V^p (p, fc; E) ^fj^'J^\, Gp^ (fc, g; K, cu) 



e/3{k,K;uj) 



(2.9) 



Since a detailed description of this formalism is already given in there is no need to repeat other equations. The 
formula to calculate the Scheerbaum factor for the s.p. spin-orbit potential by using the G-matrix solution is also 
given in We only repeat how we deal with the energy dependence of the quasi-potential Vya{p,<l', E) in the 

G-matrix equation (2.9). The total energy of the two interacting particles in the nuclear medium is not conserved. 
Since we only need the diagonal G-matrices for calculating s.p. potentials and the nuclear-matter properties in the 
lowest-order Brueckner theory, we simply use 



E" 



E" 



2/ic 



(2.10) 



both in V^aip, 9; E) and Ky^(p, k; E) in Eq. (2.9). The meaning and the adequacy of this procedure are discussed in 
|l4| by using a simple model. 



B. Effective meson-exchange potentials for fss2 



The EMEP at the quark level is most easily formulated in the momentum representation, by using the second- 
order perturbation theory with respect to the quark-baryon vertices. We employ the following qq interaction, which is 
obtained through the non-relativistic reduction of the one-boson exchange amplitudes in the parameter 7 = {m/2mud) 
(where m is the exchanged meson mass and mud is the up-down quark mass): 



47r 



U^'iqf^q. 



U''iqf,q,) = 



k^ + m? 
1 



- 1 
47r 



1 



2m^ , 2m^ , 

ua ua 



in ■ S 



m? , k^ + m? 

-rr I 



(cTi • k){cT2 • fc) - (1 - cs){m^ + k^)-{cri ■ era) 



fc2 + m' 



3g^ 

2m2 



{niudmy 



(cTi • n){cr2 • n) - (1 - Cqss)^n'^{cri ■ (T2) 



-in ■ S 



mudin 



(2.11) 



'^We use the G auss -Legendre 20-point quadrature formula to carry out the numerical integration for the partial-wave decom- 
position of Eq. (2.8). 
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Here k = Qf — Qi, a = (1 /2) (qrf + qi), and the quark- meson coupling constants are expressed in the operator form 
in the flavor space |18 19|]. For example, the product of the two different coupling-constant operators g and / are 
expressed as 



91 fi 



for 



singlet mesons 
octet mesons 



(2.12) 



where Aa(i) represents the Gell-Mann matrix for particle i. For the realistic description, the meson mixing between 
the flavor singlet and octet mesons is very important, which implies 



fr,' = fi cos 9 + fs sin 6* As , = -/i sin 9 + fs cos 6* As 



(2.13) 



instead of /i and /sAs in Eq. (2.12) for the PS mesons. Similar transformation is also applied to the V-meson coupling 
constants. The SU^ parameters of the EMEP coupling constants are therefore /i, /s, and 9. The S-meson exchange 
EMEP in Eq. ( 2.11 ) involves not only the attractive leading term, but also the momentum-dependent term and 
the LS term. The PS-meson exchange operator is the same as before, but the parameter cs is introduced only for the 
one-pion exchange, in order to reduce the very strong effect of the delta-function type contact term involved in the 
spin-spin interaction. The case cs = 1 corresponds to the full expression, while cs — corresponds to the case with 
no spin-spin contact term. The V-meson exchange potential is composed of the electric-type term, the magnetic-type 
term and the cross term. In the electric term, the central force generated by the w-meson exchange potential is usually 
most important, and it also includes the g^-type momentum-dependent term. As to the introduction of the vector- 
meson EMEP to the quark model, some discussion already addressed the problem of double counting, especially with 
the strong short-range repulsion from the u> meson [ pO[ . We will not discuss this problem here, but take a standpoint 
to avoid this double counting problem for the short-range repulsion and the LS force, by simply choosing appropriate 
coupling constants for vector mesons. The magnetic term is usually important f or th e isovector p meson, and yields 
the spin-spin, tensor and QLS terms in the standard OBEP. The choice in Eq. ( 2.11 ) is to keep only the QLS term 
with the spin-spin term proportional to L"^, the reason for which is discussed below. Finally, the cross term between 
the electric and magnetic coupling constants leads to the LS force for the qq interaction. The antisymmetric LS 
{LS'--^) force with S = {cri - (T2)/2 is not generated from EMEP, because the flavor operator in Eq. ( |2.12| ) is the 
Gell- Man n matrix and also because the mass difference between the up-down and strange quark masses is ignored in 
Eq. ( pll ). 

We should keep in mind that these EMEP are by no means a theoretical consequence of the real meson-exchange 
processes taking place between quarks. First of all, the static approximation used to derive the meson-exchange 
potentials between quarks is not permissible, since the masses of S mesons and V mesons are more than twice as 
heavy as the quark mass niud ~ 300 - 400 MeV. Since the parameter 7 is not small, the non-relativistic reduction 
is not justified. Also, the very strong S-meson central attraction is just a replacement of the real processes of the 
2tt exchange, the irp exchange, the A excitations and so forth. The V m esons are supposed to behave as composite 
particles of the (qq) pairs . Furthermore, the choice of terms in Eq. (2.11) is quite ad hoc and phenomenological. We 
should consider Eq. ( 2.11 ) as an effective interaction to simulate the residual interaction between quarks, which is not 
taken into account by the FB interaction. 



The calculation of the basic Born kernel in Eq. (2^) for each term of Eq. ( 2.11 ) becomes rather involved, if we use 
the standard technique of calculating the exchange kernel via the generator-coordinate kernel (GCM kernel). This 
becomes especially tedious, when the qq interaction involves the non-static q^ dependence and the second-order term 
of q as in the QLS force. We have developed in Q a new technique to calculate the Born kernel directly from the 
two-body interaction in the momentum representation. In this technique, there is no need to calculate the GCM 
kernel. Since the final expression is rather lengthy for the exchange kernel, it is relegated to Appendix A. Here we 
only show the direct term, which is particularly useful to see the main characteristics of the EMEP introduced in the 
present model: 



Ml{qf,q,)^g^ 



4n 



- 1 



1 



47r 



An 



fc2 



rC 3 



1 



2(3to„<j)^ 



9 

2P 



2(3m„<i) 



1 



(cTi • k){(T2 • fc) - (1 - cs){m^ + A;2)-(cri • 0-3) 

, _9_ 

^ 2fo2 



2(3m„d) 



{Sniudm) 



-n? k'^\ 3 
(cTi • n){(T2 • n) - (1 - Cqss) ( Y + ~ ) ('^1 • '^2) + TTi^ki X fe] • [era x k] 
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262' 



5 



Sruudm 



(2.14) 



Here Xqjj^ represents the spin- flavor factors related to the spin-flavor operators in Eq. (2.11). The Gaussian factor 
exp{— (6fe)^/3} appearing in Eq. ( ^.14 ) represents the form factor effect of the (Os)^ cluster wave functions. The finite 



size effect of the baryons also appears as the constant zero-point oscillation terms accompanied with the q terms, 
appearing in the S- and V-meson contributions. For the QLS force, the same effect appears as the tensor force having 
the form [cri x fc] • [(T2 x fe] . The magnitude of this term is about one third, if we compare this with the strengt h from 
the original tensor term appearing at the level of qq interaction. The advantage of using the QLS force in Eq. (|ll]), 
instead of the tensor force, is that we can avoid the Tr-p cancellation of the tensor force for the coupling term of the 
S and D waves. The ei parameter of the NN interaction is very sensitive to this coupling strength. 



C. Determination of parameters 

We have four quark-model parameters; the harmonic-oscillator width parameter b for the (Sg) clusters, the up-down 
quark mass m„(ij the strength of the quark-gluon coupling constant as, and the mass ratio of the strange to up-down 
quarks A — (ms/rriud)- A reasonable range of the values for these parameters in the present framework is 5 = 0.5 - 0.6 
fm, m„rf — 300 - 400 MeV/c^, as ^ 2, and A = 1.2 - 1.7. Note that we are dealing with the constituent quark model 
with explicit mesonic degree of freedom. The size of the system determined from the (3g) wave function with b (the 
rms radius of the {3q) system is equal to b) is related to the quark distribution, which determines the range in which 
the effect of the FB interaction plays an essential role through the quark-exchange kernel. The internal energies of 
the clusters should be calculated from the same Hamiltonian as used in the two-baryon system, and contain not only 
the quark contribution but also various EMEP contributions. The value of as is naturally correlated with 6, mud, 
and other EMEP parameters. This implies that as in our framework is merely a parameter, and has very little to do 
with the real quark-gluon coupling constant of QCD. 

For the EMEP part, we have three parameters /i, /§, and 9 for each of the S, PS, Ve (vector-electric) and Vm 
(vector-magnetic) terms. It is convenient to use the coupling constants at the baryon level, in order to compare our 
result with t he pr edictio ns b y other OBEP models. These are related to the coupling constants at the quark level 



used in Eqs. (2.11) and (2.14) through a simple relationship 



S S PS PS ^ 

A = 3gi , /s = 58 , /i = /i 7 /s = g/s , 

nVc Q fVc re /'Vm fm rVm /c\ -i r\ 



Through this replacement, the leading term for each meson in Eq. (2.14) precisely coincides with that of the OBEP 
with Gaussian form factors. In the present framework, the S-meson masses are also considered to be free parameters 
within some appropriate ranges. We further introduce three extra parameters, cs the strength factor for the delta- 
function type spin-spin contact term of the one-pion exchange potential (OPEP), Cqss the strength factor for the 
spin-spin term of the QLS force, and CqT the strength factor for the tensor term of the FB interaction. These 
parameters are introduced to improve the fit of the NN phase shifts to the empirical data, as is discussed below. 

We determine these parameters by fitting the most recent result of the phase shift analysis SP99 |^ for the np 
scattering with the partial waves J < 2 and the incident energies Tiab < 350 MeV, under the constraint of the deuteron 
binding energy and the ^5*0 NN scattering length, as well as to reproduce the available data for the low-energy YN 
total cross sections. The result is shown in Table I. The parameters of the previous model FSS are also shown for 
comparison. The value used in the parameter search is defined through 

^={^E(^r'-r'')'| , (2.16) 

where no experimental error bars are employed because the energy-dependent solution of the phase-shift analysis does 
not give them. In Eq. ( 2.16| ) the sum over i = 1 - A'^ is with respect to various angular momenta and energies, and the 
mixing parameters, ei and £2, are also included in the unit of degrees. The value a/x^ therefore gives some measure 
for the averaged deviation of the calculated phase shifts from the empirical values. Using the parameter set in Table 
I, we have obtained -^/^ = 0.655° for the np scattering. The best solution in our previous models is \/x^ ~ 3° in 
FSS. Since the present model fss2 is a renovated version of FSS, we summarize in the following only the changes and 
new points of fss2, in comparison with the model FSS: 
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1) In the original expression of the meson-exchange potentials between quarks, the momentum-dependent Bryan- 
Scott term appears in the combination of — k"^ /A for the S meson and 3q^ — for the V meson. We find 
that these terms (usually replaced by — —m?) play a rather characterless role in making the whole 
interaction slightly repulsive. With these terms, the energy-dependence of the ^Sq and ^Si phase shifts becomes 
too strong to keep the value of h in the reasonable range. (The value of h turns out to be too small, about 
h ^ 0.4 fm to compensate the strong energy dependence.) We therefore drop all these A;^/4 terms in the present 
calculation. 

2) We ignore the QLS force from the S-mesons, since it is very weak. The S-meson EMEP therefore consists of 



the leading term with —1 in Eq. (2.14), the momentum-dependent Bryan-Scott term and the LS term. This LS 
term yields an appreciable contribution at medium and higher energies, which consequently reduces the value 
of h from the previous value > 0.6 fm to a smaller value ~ 0.56 fm. 

3) The reduction of the spin-spin contact term for the PS mesons is introduced only for the pion with the smallest 
mass. For the other heavier PS mesons, we assume the full strength factor q = 1. The reduction from 1 for the 
pion improves the fit of the NN ^Pi phase shift to a great extent. (Otherwise, the repulsion at higher energies 
is insufficient for this partial wave.) We introduce cs only for pion, since the effect of the present (3(7)-cluster 
folding corresponds to a very low value of the cut-off mass A ~ 800 - 900 MeV for the pion form factor in OBEP. 
It is well known that such a low value of A converts even the sign of the medium-range part of the OPEP if 
the full strength of the contact term is introduced. The factor ca < 1 also reduces the very strong repulsion 
generated from the one pion spin-spin contact term for the S'-wave states of the NN system. In the present 
framework, this repulsion is almost 300 MeV if cs = 1 is assumed. Furthermore, the value of cs has a strong 
influence on the internal energies of single baryons. It reduces the very large contribution of the pion to the 
A^-A and A-S mass difference, the latter helping us to keep A = (ms/mud) at the moderate value. (Otherwise, 
we obtain A ~ 1.) If we do not introduce cs and the parameters Cqss, CqT discussed below, the \/x^ value cannot 
be improved by more than 1.5°. The contribution of t] and 77' mesons was necessary in the previous models in 
order to make the central force relatively more repulsive than the ^S* central force. In the present framework, 
it turns out that the introduction of these rj mesons is not convenient for the subtle balance of the central and 
tensor forces. We therefore take out all these 77-meson contributions. The well-known too strong repulsion of the 
A''A^ central force from the color-magnetic interaction of the FB interaction |^l| is remedied by assuming two 
different masses for the isovector 5 meson, i.e., ms — 720 MeV/c^ for the NN system and ms — 846 MeV/c^ 
for the YN system (see comment 3) in Table I). 

4) As is discussed at the end of the preceding subsection, the present model fss2 is the QLS dominant model. 
This implies that we use the QLS force to reduce the too strong OPEP tensor force, instead of the tensor force 
itself. The main reason for this choice is that the NN mixing parameter ei is very difRcult to reproduce if the 
cancellation of the one pion tensor force and the p-meson tensor force is too strong for the S'-wave and ZJ-wave 
coupling. Another question is how this QLS force is incorporated into the model. We find that the QLS spin- 
spin term n^(cri ■ (T2) in Eq. ( B.ll ) plays a favorable role in improving the fit of the NN phase shifts. This term 



corresponds to the (cri • cr2)i" term in the Hamada- Johnstone potential Since the full introduction of this 
term results in too vigorous behavior, we introduce a reduction factor Cqss^ which turns out around Cqss 0.6. 
The two-pole formula for the p-meson exchange potential, introduced in p^ , is found to give a favorable result. 
We further find that the short-range tensor force is still too weak. We avoid this difficulty simply by increasing 
the strength of the tensor term of the FB interaction with the factor Cqx- The value Cqx ^ 3 seems to be 
reasonable. If we carry out the parameter search with CqT = 1, the value of a/x^ cannot be improved by more 
than 1.3° ^ 1.0°, mainly due to the disagreement of ei. We should note, however, that the introduction of the 
V mesons is a rather minor change from our previous models. With the exception of f^™ — 2.577, the V-meson 
coupling constants in Table I are around one, which is less than half of the coupling constants in the standard 
OBEP. In particular, the isospin dependent LS force from the p meson is exactly zero, since /g'^ is fixed at zero. 
The short-range repulsion in the NN interaction is still mainly described by the color-magnetic term of the FB 
interaction. The dominant effect of the V mesons is almost solely the p meson QLS force, which is the reason 
we call fss2 the QLS dominant model. 

5) The following five parameters in Table I are directly related to the reproduction of the low-energy YN cross 
sections; A = {ms/mud), d^, df, w^, and ttIk,. Among them, the angle of the singlet-octet meson mixing 9^ of 
the S mesons are used to control the relative strength of the central attraction of the NN and YN interactions. 
It was found before [Q that, once the 9^ is determined to fit the low-energy Ap cross section data, the attraction 
of the T,N{I — 3/2) channel is too strong and the S+p total cross sections are overestimated. We therefore use a 
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larger value for 6^ (which is denoted by 6'f ) only for the T,N{I = 3/2) channel in order to reduce the attraction, 
which is the same prescription employed in the previous models 



6) The largest ambiguity for determining the parameters related to the YN interaction lies in the strength of 
the central attraction in the I]N{I = 1/2) "^^i channel 0. If the phase-shift rise of the '^Si state is less than 
30°, the low-energy elastic total cross section becomes too small. If this attraction is too strong, as in 
RGM-F the "^Si phase shift shows a sudden decrease from 180° to 60 - 90°, and the behavior of the Ap 
total cross sections at the EiV threshold becomes a round peak, instead of the cusp structure |Q . Furthermore, 
the strength of the central attraction plays a crucial role even for the odd-parity state. The T,N{I — 1/2) 
•^Pi phase shift is attractive due to the exchange kinetic-energy kernel; i.e., the effect of the Pauli principle 
[ p^ . This attraction is reinforced by the LS force in the diagonal channel, and also by the LS^^'' force acting 
between this channel and the ^Pi channel. This channel coupling also takes place between the 'SN{I — 1/2) 
channel and the AA^ channel. This channel coupling is mainly determined by the strength of the LS^~^ force, 
which is directly related to the magnitude of as, but also considerably influenced by the strength of the central 
attraction in the TtN{I = 1/2) channel. In [Q, we have clarified that the central attraction of the previous 
models RGM-F and FSS is so strong that the 'SN{I — 1/2) ^Pi resonance is moved to the AA^ ^Pi channel. 
The consequence of this behavior is the strong enhancement of the Ap total cross sections in the cusp region. 
On the contrary, the P-wave coupling in the model RGM-H is less strong, and the agreement of the Ap total 
cross sections is much better. (See Fig. 10(a) in [|j and Table II in Q.) Here we assume that the resonance stays 
in the original T,N{I = 1/2) ^Pi channel, and try to find the parameter set which gives the maximum strength 
of the EN (I = 1/2) central attraction. In practice, we assume y/2/TTasx^mudc'^ = 440 MeV {x = {fi/mudcb) is 
the ratio of the Compton wave-length of the up-down quarks to 6) as in RGM-F and FSS,0 and adjust the value 
of ms for the YN interaction, independently of the value in the case of A^A^ interaction. If we use a smaller 
value for m^, the T,N{I = 3/2) ^Sq state becomes more attractive and the T,N{I — 1/2) '^S'l state becomes less 
attractive. 

7) Another important change from the previous models FSS and RGM-H is the relative strength of the ^5*0 and 
•^5*1 attraction in the AA^ interaction. The maximum phase-shift values of the ^So and ^S'l states in these models 
are about 46° and 16°, respectively, around pA ^ 200 AIcV/c. The big difference of almost 30° is known to be 
unfavorable for the description of the s-shell A-hypernuclei. Detailed few-body calculations for these hypernuclei 
have recently been carried out by several groups | 
effective AA^ interactions, the effect of the Y,N channel coupling is usually renormalized. These calculations 
imply that the phase-shift difference of a little less than 10° seems to be most appropriate. We follow this 
suggestion and adjust the strength of the AA^ attraction such that the ^5*0 and ^5*1 phase-shift difference is less 
than 10° and the low-energy Ap cross sections are correctly reproduced. We can use the K-meson mass to adjust 
this phase-shift difference. Namely, if is smaller, then the AA^ ^Sa phase shift becomes more attractive and 
the ^Si phase shift becomes less attractive. 

In order to give an outline of the framework, we summarize the difference of FSS and fss2 in Table II, with respect 
to the meson species and interaction types of EMEP included in the models. Table HI shows the quark and EMEP 
contributions to the baryon mass difference between A^ and A {AEn-a — -E^* — £^]^*), and the mass difference between 
A and E (AE'a-s — E^^ — £'X'*), calculated in the isospin basis. We note that various meson contributions largely 
cancel each other and the net contribution is roughly given by the quark contribution from the color-magnetic term 
of the FB interaction. 



D. Calculation in the particle basis 

In this subsection we discuss some new features required in the calculation in the particle basis. Three different 
types of calculations are carried out in this paper. 

1) calculation in the isospin basis 



^This value corresponds to assuming the N-A mass difference 293.3 MeV only by the FB interaction, as seen from Table III. 
If we use the as value about 1.3 times larger, the transition of the P-wave resonance to the A7V ^Pi channel takes place in the 
present model. 



25-p8| by using various effective AN interactions. In these 
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2) calculation in the particle basis without the Coulomb force 

3) calculation in the particle basis with the Coulomb force 

For the NN interaction, the calculation in the particle basis is rather straightforward. We use the empirical baryon 
masses listed in Table |V| and evaluate spin-flavor factors for the charged pion and the neutral pion separately in 
the isospin representation. The other spin-flavor factors for heavier mesons and the FB interaction are generated in 
the simple isospin relations. The Coulomb force is introduced at the quark level by using the quark charges. The 
exchange Coulomb kernel has the same structure as the color- Coulombic term of the FB interaction. 

Only complexity arises when we solve the LS-RGM equation in the momentum representation. The standard 
technique by Vincent and Phatak p9[ is employed to solve the Lippmann-Schwinger equation in the momentum 
representation, including the Coulomb force. This technique requires introducing a cut-off radius Rc for the Coulomb 
interaction. In the RGM formalism, we have to introduce this cut-off at the quark level, in order to avoid violating 
the Pauli principle. The two-body Coulomb force assumed in the present calculation is therefore written as 

U!^"" = Q^QJe'—Q{Rc - nj) , (2.17) 

rij 

where Q is the Heaviside step function and Qi, Qj — 2/3 for the up quark and —1/3 for the down and strange quarks. 
The Coulomb contribution to the internal energies becomes zero for the proton and S+. More explicitly, this can be 
given by 



ES^t = X^j^^axmudc" (^1 - e iC^'' ) j , (2.18) 

where a = {e^ /fic) ~ 1/137 is the hyperfine coupling constant and the direct spin- flavor factor is expressed as 
= X]i=i 2 ~ l/3)/2 — 1/3] in terms of the total charge Zi of the i-th baryon. The basic Born kernel for 

the direct Coulomb term reads 

M^^Qf^q,) = Z^Z^eH^Rc^ (^-^ sin ^ j e-^('"=)' with fc = |q/ - q^\ , (2.19) 

which corresponds to the direct Coulomb potential 



VD{r) = Z^Z^e^-^ jerf (V70 " \ [erf (^/^(r -f Rc)) + erf (^^(r ~ i?c))]| 



(2.20) 



Here erf (a;) = i^j^) e~^^dt stands for the error function and ^ = jiv = (3/46^). The exchange Coulomb 
kernel is also slightly modified from the exact Coulomb kernel. This is given in Appendix A, together with other 
EMEP kernels. The value Rc should be sufficiently large to be free from any nuclear effect beyond Rc- Then the 
final ^-matrix is calculated from the condition that the wave function obtained by solving the Lippmann-Schwinger 
equation with the modified Coulomb force is smoothly connected to the asymptotic Coulomb wave function. We take 
Rc = 9 fm, although a much smaller value seems to be sufficient. Note that, even in the n-p and nn systems, we have 
small contributions from the Coulomb interaction through the exchange Coulomb kernel. The difference between the 
calculations 2) and 3) for the system of chargeless particles implies this effect. 

For the YN interaction, more consideration is required for the treatment of the threshold energies. We note that 
the mass difference of S~ and E"*" is about 8 MeV and is fairly large. Figure |l| shows the comparison of threshold 
relations in the isospin and particle bases, evaluated in the non-relativistic kinematics. The Ap system has the total 
charge Q = \ and the system (5 = 0. The direct Coulomb term exists only in the S+p channel and the S~p 
channel. The EMEP contribution to the A - E mass difference is given in Table V, both in the isospin basis and in 
the particle basis. We find rather large cancellation between the neutral and charged pion contributions to the Yi^ 
- TP mass difference. The pion-Coulomb correction yields the calculated threshold energies given in Table VI in the 
non-relativistic kinematics. Apparently the empirical mass difference is not precisely reproduced.^ When we use the 
relativistic kinematics, the empirical threshold energies are defined by E'/^ = (p™)^/2/iinc (Mine is the non-relativistic 
reduced mass of the incident channel), which are used in the present non-relativistic model in the cm. system. Table 



^The up-down quark-mass difference does not help, because the total hypercharge is conserved in the two-baryon systems. 
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VII shows these energies in the column of relativistic E^^^. The difference between these E^^^ and the calculated 
threshold energies E^^^ is given in the last column. For the S~p - - An system, E^^^ are calculated from those in 
the An - E°n - E~p system. The disagreement of the threshold energies between the calculation and the experiment 
is a common feature of any microscopic model. Fortunately, we have a nice method to remedy this flaw without 
violating the Pauli principle. As discussed in ||l^ in detail, we only need to add a small correction term AG to the 
exchange kernel, in order to use the empirical threshold energies E^^^ listed in Table VII. The same technique is also 
applied to the reduced mass corrections. 



III. RESULTS AND DISCUSSIONS 



A. The NN result 

Figures ||(a) - ||(i) compare the np phase shifts and the mixing angles ej predicted by fss2 with the recent phase-shift 
analysis SP99 by Arndt et al. Q The parameter search and the calculation of phase-shift parameters in this subsection 
are all carried out in the isospin basis. For comparison, the previous results by FSS are also shown with the dotted 
curves. Here we examine the partial waves up to J = 4 in the energy range Tjab = - 800 MeV. For energies higher 
than 300 MeV, the inelasticity parameters of SP99 are given for a measure of possible deviations of the phase-shift 
values in the single-channel calculation. The phase shift is greatly improved by the QLS component. Even in 
the other partial waves, the improvement of the phase-shift parameters is usually achieved. This includes 1) ^Pq, ^Pi, 
and ^G4 phase shifts, 2) ^Si, ^So, ^Pi, ^Fs, and phase shifts at higher energies Tiab = 400 - 800 MeV, and 3) 
some improvement in phase shift and £2 mixing parameter. On the other hand, ^P2 and ^D^ phase shifts turn out 
worse and phase shift is not much improved. The disagreement of the ^D^ phase shift and the deviation of the 
"^Di phase shift at the higher energies imply that our description of the central, tensor and LS forces in the ^E states 
requires further improvement. The insufhciency in the partial waves is probably related to the imbalance of the 
central force and the LS force in the short-range region. The decomposition of the "^Pj phase shifts to the central, LS 
and tensor components, shown in Fig. ||, implies that the central force is too repulsive at higher energies Tiab > 400 
- 500 MeV. It should be noted that whenever the discrepancy of the phase-shift parameters between the calculation 
and the experiment is large, the inelasticity parameters are also very large. In particular, the inelasticity parameters 
of the ^P2, ^D2, and ^F^ states rise very rapidly as the energy increases, and reach more than 20° at Tiab = 800 
MeV. The elastic phase shift for each of these states shows a dispersion-like resonance behavior at the energy range 
from 500 MeV to 800 MeV. These are the well-know di-baryon resonances directly related to the AA^ threshold in 
the isos pin / — 1 channel. The present single-channel calculation is not capable of describing these resonances. 



Table VIII tabulates the values of phase-shift parameters in the energy range Tiab = 25 - 300 MeV, in comparison 
with those of SP9 9 [p| and other meson-exchange models. The results of OBEP, Paris and Bonn potentials are cited 
from Table 5.2 in ||31[| . The partial waves only up to J = 2 are considered. If we calculate the values using these 
numbers, we obtain \px^— 0.59, 1.10, 1.40 and 1.32 for fss2, OBEP, Paris and Bonn, respectively. The reason we 
get such results is as follows. In the meson exchange models, the accuracy of the low energy phase shifts is less 
than 0.2°, and the agreement with the experiment is excellent. However, in higher energies the deviation from the 
experiment increases, and in some particular partial waves like ^5*0 and ■^Pq states, it becomes more than 2°. In the 
Paris potential, the ^6*0 phase shift is apparently too repulsive. This is, however, because the parameters of the Paris 
potential is determined by the fit to the phase shifts, and the correction due to the CSB is probably not taken into 
account in the numbers given here. Every model has its own weak points. For example, the tensor force of the Bonn 
potential is usually very weak, which is reflected in the e\ parameter and in the too attractive behavior of the '^Pq 
phase shift. (However, the recent CD Bonn potential ||T^ fits the A^A^ phase-shift parameters in the non-relativistic 
energies almost perfectly, with various possible corrections taken into account.) The weak point of our model lies in 
the ^P2 and ^D-^ phase shifts at the intermediate and higher energies Tiab = 300 - 800 MeV. The empirical ^P2 phase 
shift gradually decreases if we ignore the weak dispersion-like behavior. Our result, however, decreases too rapidly. 
Our ^1)3 phase shift is too attractive by 4° - 6°. 

Figures]^ and |^ illustrate the fss2 predictions of the differential cross sections {da/d£L in mb/sr) and the polarizations 
{P{0)) for the elastic np scattering, in comparison with experiment The same observables for the elastic pp 
scattering are plotted in Figs. ^ and 0. The final calculation of these observables in this paper is carried out in the 
particle basis with the full Coulomb force incorporated. The corresponding figures by our previous model FSS are 
given in Figs. 1, 2 of and Figs. 2, 3 of 0. (Note that the plot of the differential cross sections at higher energies 
in j^] is given in the log scale.) We find some improvement in the differential cross sections. First, the previous 
overestimation of the np differential cross sections at the forward angle at Tiab = 320 MeV is corrected. Secondly, 
the bump structure of the np differential cross sections around 0cm = 130° at energies Tiab = 300 - 800 MeV has 
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disappeared. The overestimation of the pp differential cross sections at 6'cm = 10° - 30° at energies Tiab = 140 - 400 
MeV is improved. However, the essential difficulties of FSS and RGM-H, namely the oscillatory behavior of the np 
polarization around Ocm ^ 110° and that of the pp polarization around the symmetric angle 6'cm = 90° for higher 
energies Tiab > 400 MeV are not resolved. Furthermore, the pp differential cross sections show a deep dip at angles 
6'cm < 30° and > 150° for Tiab > 500 MeV. The low-energy pp cross sections at 6'cm = 90° for Tiab < 100 MeV are 
still overestimated. 

In order to find a possible reason for the unfavorable oscillations of our polarizations, we show in Fig. || the five 
independent pp invariant amplitudes at the highest energy Tiab = 800 MeV. They are composed of the real and 
imaginary parts of 170 (spin-independent central), (LS), ft.„ ( {cri ■ n)((T2 • n)-type tensor), /ij. ( {cri ■ k)[(j2 ■ A;)-type 
tensor), and hp ( [cti ■ P){(T2 ■ P)-type tensor) invariant amplitudes. In Fig. |^ the Coulomb force is neglected in 
the predictions by the Paris potential. The result by SP99 is calculated using only the real parts of the empirical 
phase-shift parameters. If we recall that the polarization is given by the cross term contribution of the central, LS, 
and tensor invariant amplitudes (i.e., P{9) = 29rn [{go + hn){ho)*], see Eq. (2.32) of 0) we find that the disagreement 
in 5m hn and SRe ho is most serious. Since the oscillatory behavior of 9m /i„ in SP99 also appears in 5m hk and 
5m hp, it is possible that this is an oscillation caused by the NN - AN channel coupling through the one pion 
spin-spin and tensor forces. Figure || also shows the reason for the underestimation of the differential cross sections 
at 6'cm < 30°. Namely, The imaginary part of go is too small both for fss2 and the Paris potential, and the real part 
of go is strongly reduced in fss2. 

Another application of the invariant amplitudes is the p prescription for calculating the s.p. potentials of the 
nucleons and hyperons in nuclear matter. It is discussed in that the s.p. potentials predicted by the model FSS 
in the G-matrix calculation show fairly strong attractive behavior in the momentum interval qi =5-20 fm~^ for all 
the baryons. In particular, UN{qi) in the continuous prescription becomes almost —80 MeV at qi = 10 fm"^. This 
momentum interval corresponds to the incident energy range Tiab = 500 MeV -8 GeV in the NN scattering. The 
prescription is a convenient way to evaluate the s.p. potentials in the asymptotic momentum region in terms of 
the spin- independent invariant amplitude at the forward angle go{d — 0). Since the present model fss2 incorporates 
the momentum-dependent Bryan-Scott term, the asymptotic behavior of the s.p. potentials in the large momentum 
region is improved. We can see this in Fig. ||, where the s.p. potentials of N , A, and E calculated in the G-matrix 
approach are shown in the momentum range gi = - 10 fm~^. Figures |9|(a) and (b) show the result in the QTQ 
prescription, and Figs.||(c) and (d) in the continuous choice for intermediate spectra.^ Figures ||(a) and (c) show the 
real part of Usiqi), and Figs.^(b) and (d) the imaginary part. In Figs.^(c) and (d), the sohd curves for the nucleon 
s.p. potential are compared with the results by the p prescription with respect to the T-matrices of fss2, the Paris 
potential |5^, and the empirical phase shifts SP99 ||]. The partial waves up to J < 8 are included in fss2 and the 
Paris potential, and J < 7 in SP99. The momentum points calculated correspond to the energies Tab — 100, 200, 400, 
800, and 1,600 MeV. We find that the real part of UN{qi) nicely reproduces the result of the G-matrix calculation 
even at such a low energy as Tiab = 100 MeV. On the other hand, the imaginary part by the t°*p prescription usually 
overestimates the exact result especially at the lower energies. 

B. Deuteron properties and effective range parameters 

The deuteron properties are calculated by solving the LS-RGM equation with respect to the relative wave functions 
/o(fc) and /2(fc) in the momentum representation (see Appendix C). The properly normalized wave functions in the 
Schodinger picture are not fi{k) but Fi — \/N ft, where N represents the normalization kernel [Q. The iS-wave 
and T'-wave wave functions in the coordinate representation, u{R) and w{R), are then obtained from the inverse 
Fourier transform of Fi(k). This process is most easily carried out by expanding Fi{k) in a series of Yukawa functions 
^J^fnklik'^ -|- 7j^) in the momentum representation (see Appendix D in |l|]). We choose 7^ = 7 + {j — 1)70 with 
7o = 0.9 fm~ and j = 1 - 11. The 7 is the S'-matrix pole q — —17, from which the deuteron energy is most 
accurately calculated by using the relativistic relation 

M„ + Mp - Ed = V^M„2-^2 + ^Mp2_-y2 . (3.1) 



*In Fig.g we have employed the approximate angular integration given in Eq. (A. 8) of [g], while in Fig. ^ the numerical 
integration over the angle Q2 in Eq. (A. 6) of is explicitly carried out. 
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Figure [T^ shows the deuteron wave functions of fss2 in the coordinate and momentum representations, compared 
with those of the Bonn model-C potential (dotted curves)0. We find that the difference between the two models 
is very small. Table IX compares various deuteron properties calculated in three different schemes. They are also 
compared with the empirical values and the predictions by the Bonn model-C potential. The final value of the 
deuteron binding energy for fss2 is — 2.2309 MeV. If we use the non-rclativistic energy expression,^ Cd = (7^/Mjv) 
for 7^ = 0.05376157 fm"^ in the fuU calculation, we obtain e<j = 2.2295 MeV and the difference is 1.4 keV. The 
differences within the deuteron parameters calculated in the three different schemes are very small, except for the 
binding energy e^- In particular, the exchange Coulomb kernel due to the exact antisymmetrization at the quark level 
gives an attractive effect to the binding energy, and increases by 4.8 keV. This is even larger than the relativistic 
correction included in Eq. (3.1). The deuteron D-state probability is Pd = 5.49 % in fss2, which is slightly smaller 
than 5.88 % in FSS [Q. These values are rather close to the value Pjj = 5.60 % obtained by the Bonn model-C 
potential pl[ |. The asymptotic D/S state ratio 77 and the rms radius are very well reproduced. On the other hand, 
the quadrupole moment is too small by about 5 - 6%. There are some calculations [^3|,^ which claim that the effect of 
the meson-exchange currents on the dueteron quadrupole moment is as large as AQd = 0.01 fm^. It is noteworthy that 
the Bonn model-C almost reproduces the correct quadrupole moment, in spite of the fact that the I?-state probability 
is very close to ours. (On the other hand, the quadrupole moment of CD-Bonn [|ll| is Qd — 0.270 fm^ with a smaller 
value Pd — 4.85%.) For the magnetic moment, precise comparison with the experimental value requires a careful 
estimation of various corrections arising from the meson-exchange currents and the relativistic effect, etc. 

Table ^ lists the S- and P-wave effective range parameters for the NN system, calculated in the three schemes. 
Since the pion-Coulomb correction is not sufficient to explain the full CIB effect existing in the np and pp ^Sg states, 
a simple prescription to multiply the flavor-singlet S'-meson coupling constant ff by a factor 0.9949 is adopted to 
reduce the too large attraction of the pp central force. (This prescription is applied only to the calculation in the 
particle basis.) The underlined values of a in Table ^ indicate that they are fitted to the experimental values. We find 
that the pion-Coulomb correction in the np ^Sq state has a rather large effect on the scattering length parameter a. 
The value a = —23.76 fm in the isospin basis changes to a = —27.38 fm due to the effect of the pion mass correction 
and the explicit use of the neutron and proton masses. It further changes to a = —27.87 fm due to the small effect 
of the exchange Coulomb kerne l. T hese changes however should be carefully reexamined by readjusting the binding 
energy of the deuteron in Table IX , Wc did not carry out this program, since the reduction of ff to fit these values to 
the empirical value a = —23.748 ± 0.010 fm does not help much to reproduce the CIB of the pp channel anyway. We 
have to say that the improvement of the NN 5*- wave effective range parameters in the particle basis calculation is not 
excellent, in spite of the large effort expended in incorporating the pion-Coulomb correction in the microscopic RGM 
formalism. This shortcoming might be related to the insufficient description of the low-energy pp differential cross 
sections around 6'cm ^ 90°, observed in Fig. ^. It was also pointed out by the Nijmegen group that the Coulomb 
phase shift should be improved by the effects of two-photon exchange, vacuum polarization and magnetic moment 
interactions, in order to describe the ^5*0 phase shift precisely at energies less than 30 MeV. These effects are not 
incorporated in the present calculation. The P-wave effective range parameters are also given in Table in order to 
compare with a number of empirical predictions. The parameters of ^P2 state are not given, since the effective range 
expansion of this partial wave requires a correction term related to the accidental p^ low-energy behavior of OPEP 

m. 



C. S p system 

Figure |l^ displays the S- and P-wave phase shifts of the T,N{I — 3/2) system, calculated in the isospin basis. The 
results given by FSS (dashed curves) and RGM-H (dotted curves) are also shown for comparison. The and '^O 
states of the T,N{I = 3/2) system belong to the (22) component in the flavor SU3 representation, which is common 
with the and '^O states of the NN system |l^. The phase-shift behavior of these partial waves therefore resembles 
that of the NN system, as long as the effect of the SU3 symmetry breaking is not significant. Figure ^ shows that 
the attraction in the ^Sq state is much weaker than that of the NN system, and the phase-shift peak is about 26° 
around ps — 200 MeV/c. The "^Pj phase shifts show the characteristic energy dependence observed in the NN phase 
shifts, which is caused by competition among the central, tensor and LS forces. The appreciable difference in the 



^The results of the Bonn model-C potential in Fig. ^ and in Table ^ are based on the parameterized deuteron wave functions 
given in Table C.4 of Q. 
®In Table the value of in the isospin basis is calculated using this non-relativistic formula. 
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three models fss2, FSS and RGM-H appears only in the state. It is discussed in |4|] that the attractive behavior of 
the phase shift is closely related to the magnitude of the S+p polarization P{9) at the intermediate energies. The 
more attractive, the larger P{9 = 90°). Since the attraction of fss2 is just between FSS and RGM-F (see Fig.jlWb)), 
the polarization curve at = 450 MeV/c falls between the two curves given by FSS and RGM-H (see Fig. 4 of |142[|). 
This implies our quark-model prediction P{9 = 90°) = 0.1 - 0.2 at ps — 450 MeV/c. 

On the other hand, the '^E and states of the T,N{I = 3/2) system have the (30) symmetry of the flavor SU^ 
We have no information on the properties of these states from the NN interaction. In our present framework of the 
quark model, there is very little ambiguity in the phase shifts of these states, as can be seen from Fig. |ll|. Since the 
configuration (Os)^ in the ^S'l state is almost forbidden by the effect of the Pauli principle, the interaction in the ^S'l 
state is strongly repulsive. This property is common in all of our models, and the phase shifts predicted by fss2, FSS 
and RGM-H are almost the same. On the other hand, the ^Pi phase shift is weakly attractive, which is caused by 
the exchange kinetic-energy kernel due to the Pauli principle. This property is also common to all three models. The 
Nijmegen hard core models D and F give a resonance in the phase shift of this channel, which induces enhancement 
in the E+p differential cross sections at the forward and backward angles. (See Fig. 6 in ||^.) Such behavior is not 
found in any of our quark models. 

The effective range parameters of the YN scattering in the single-channel analysis are given in Table ^ with some 
empirical values. For the E+p system, the empirical values given in |ll3|] should be compared with the results in the 
particle-basis calculation including the Coulomb force. We find a reasonable agreement both in the ^Sq and ^5*1 states. 



D. AN system 

The total cross section for the Ap elastic scattering predicted by fss2 in the isospin basis is displayed in Fig. ^(a) , 
together with the previous result given by FSS. The new model and FSS reproduce experimental data equally well 
in the low momentum region pA < 300 MeV/c. The total cross section has a cusp structure at the 'EN{I = 1/2) 
threshold, which is due to the strong AA^ - I]N{I = 1/2) ^Si - '^Di coupling caused by the OPEP tensor force. To 
see the dominant role of the tensor force of the 'SN - AN coupling in more detail, we show the '^Si and ^Sq phase 
shifts in Fig. ^ which are calculated by fss2 (left) and FSS (right) in the isospin basis. A cusp structure at the 
T,N threshold is apparent in the ^Si channel, while very small in the ^Sq channel. In both models, attraction in the 
^5*0 state is stronger than that in the "^Si state. However, the differences between the strengths of attraction in the 
^So and ^Si states vary. The old model FSS gives SC-Sq) - S{^Si) - 30° at ~ 200 MeV/c, while fss2 predicts 
6CSo) - 6{^Si) ~ 10°. The difference of S{^So) - S{^Si) - 10° is required from the few-body calculations IfUl of 
s-shell A-hypernuclei, as discussed in 7) of Sec. II C. The parameter search of fss2 is carried out under this constraint. 

As seen in Fig.^(a), FSS predicts especially large enhancement of the total cross sections around the cusp at the 
YiN threshold. This is due to a rapid increase of the AA^ '^Pi - AN ^Pi transition around the threshold. Figure 
H shows the S'-matrix S^j = r],je'^'^'^ for the AA^ - T,N{I = 1/2) ^Pi - ^Pi channel coupling for fss2 (upper), FSS 
(middle), and RGM-H (lower). The result by FSS shows that the transmission coefficient 7721, which corresponds to 
the AA^ ^Pi — > AA^ ^Pi transition, increases very rapidly around the 'SN threshold as energy increases. This increase 
of the ?72i (and the resultant decrease of the reflection coefficient 7711) is a common feature of all our models. The 
strength of the transition, however, has some model dependence. The transition is stronger in FSS than in fss2 and 
RGM-H. In particular, the resemblance of ^-matrix in fss2 and RGM-H is very outstanding. 

The behavior of the diagonal phase shifts is largely affected by the strength of the ^Pi - "^Pi channel coupling (which 
is directly reflected in 7721)- The new model fss2 and RGM-H yield a broad resonance in the SA^ '^Pi channel. On the 
other hand, FSS (and also RGM-F) with stronger channel coupling provides no resonance in this channel. Instead, 



a step-like resonance appears in the AA^ ^Pi channel. This situation is summarized in Table Xl. The location of 
the resonance is determined by the strength of the LS''"-' force and the strength of the attractive central force in 
the EN{I — 1/2) channel. In the quark model, the central attraction by the S mesons is enhanced by the exchange 
kinetic-energy kernel, which is attractive in the SA^ ^Pi channel due to the effect of the Pauli principle. The LS 
and LS"*^"-* forces also contribute to increase this attraction. In FSS and RGM-F, a single channel calculation for the 
'EN{I = 1/2) system yields a resonance in the ^Pi state. When the channel coupling to the AA^ system is introduced, 
this resonance transfers to the AA^ channel due to the strong LS^''' force generated from the FB interaction. On 
the other hand, the central force of RGM-H in the SAr(/ = 1/2) channel is weaker than in FSS (see l^siV(i/2) (^'^) 

Table II of §). The LS'^-^ force of RGM-H is also weak {y'^asx^mudc^ = 296 MeV). Accordingly, the resonance 
stays in the SA^ ^Pi channel, even if the channel coupling is incorporated. In the new model fss2, the LS term of the 
S-meson exchange is included. Since EMEP yields no LS"*^"-* force in the present framework (see Sec. II B), the role of 
the FB LS"^"^ force becomes less significant, in comparison with the very strong effect in FSS. This is the reason the 
resonance remains in the Y.N ■^Pi channel in fss2. 
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In spite of these quantitative differences in the couphng strength, the essential mechanism of the AA'' - 'EN{I =1/2) 
^Pi - '^Pi channel coupling is the same for all of our models. It is induced by the strong FB LS'*^"-' force, which directly 
connects the two SU3 configurations with the (ll)s and (11)q representations in the P-wave / = 1/2 channel. In order 
to determine the detailed phase-shift behavior of each channel, including the position of the P-wave resonance, one has 
to know the strength of the LS^^'^ force and the strength of the attractive central force in the EA^(/ = 1/2) channel. 
This is possible only by the careful analysis of rich experimental data for the Ap and S^p scattering observables in 
the SiV threshold region. 

Let us discuss another important feature of the AA^ phase shifts induced by the P-wave coupling due to the LS^^^ 
force. The AN ^Pi and ^Pi phase shifts in Fig. ^ show weakly attractive behavior over the energies below the EiV 
threshold. This is due to the dispersion-like (or step-like for the AN ^Pi state in FSS) resonance behavior with a large 
width. This attraction in the P state can be observed by the forward to backward ratio {F/ B) of the Ap differential 
cross sections, as is seen in Fig. |l|(b). Our new model fss2 and FSS give F/B > 1 below the YiN threshold, which 
implies that the P-state AA^ interaction is weakly attractive as suggested by Dalitz et al. |4^. In our models, this 
attraction originates from the strong A7V - 'EN (I = 1/2) ^Pi - '^Pi coupling due to the FB PS'^^-' force ||]. 

Some comments are in order with respect to particle-basis calculations of the Ap and An scatterings. From the 
energy spectrum of the - \He isodoublet A-hypernuclei, it is inferred that the Ap interaction is more attractive 
than the An interaction. This CSB may have its origin in the different threshold energies of the EN particle channels 
as in Fig. |l| and the Coulomb attraction in the E~p channel for the An system. The former effect increases the 
cross sections of the Ap interaction and the latter the An interaction^ However, the full calculation including the 



pion-Coulomb correction using the correct threshold energies in Table VII yields very small difference in the S-wave 
phase shifts. In the energy region up to pa = 200 MeV/c, the Ap phase shift is more attractive than the An phase 
shift only by less than 0.2° in the ^Sq state, and by less than 0.4° in the ^Si state. This can also be seen from the 
Ap and An effective range parameters tabulated in Table obtained in the three different calculational schemes. On 
the other hand, a rather large effect of CSB is found in the E~p channel as discussed in the next subsection. This is 
because the CSB effect is enhanced by the strong AA^ - I]N{I = 1/2) channel-coupling effect in the '^S'l - ^Di state. 
Figure |l5|(a) displays the enlarged picture of total AA^ cross sections in the EN threshold region, calculated in the full 
scheme 3) with the Coulomb force. The effect of different threshold energies in the particle basis is clearly observed. 
Figure |l5|(b) demonstrates the behavior of the subthreshold reaction cross sections, which shows a very strong channel 
dependence related to the small difference of the threshold energies. In contrast to this, the Coulomb effect in the 
E^p channel is found to be very small as long as the reactions from the AA^ incident channel are concerned. 

E. The S^p system 

Since the E^p system is expressed as \E^p) = —y^2/3\EN{I — 1/2)) -I- ^/T/3\EN{I = 3/2)) in the isospin basis, 
it is important to know first the phase-shift behavior of the EN {I — 1/2) system. The SU3 decomposition of the 
EN{I ^ 1/2) state 

EAr(/ = 1/2) = ^ [3(11),. - (22)] for and states, 
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EAr(/ = 1/2) = ^ [(ll)a + (03)] for and states, (3.2) 

v2 

is very useful to know the quark- model prediction for the phase-shift behavior in the isospin basis p^ . The most 
compact (Os)^ configuration in the (ll)s SU^ state is completely Pauli forbidden. This implies that the EN{I = 1/2) 
^5*0 phase shift is very repulsive due to the exchange kinetic-energy kernel On the other hand, EN{I = 1/2) 
"^Si phase shift is expected to have attraction similar to the AA^ ^5*1 phase shift, as long as the effect of the flavor 
symmetry breaking is not important. (Note that AA^ = [—(11)^ -I- (03)] /\/2 for ■^P and states.) Unfortunately, 
the last condition is applicable only to the central force, since the one pion tensor force introduces quite a few 
complexities in the channel dependence in the EN{I = 1/2) - AA^ '^S'l - '^Di channel-coupling problem. The strength 
of the EN{I = 1/2) central attraction, discussed in the preceding subsection, should therefore be examined carefully 
after this S - D wave channel coupling is properly treated. 

Figure displays the ^5*0 phase shift and the S'-matrix of the AN-EN{I = 1/2) - '^Di coupled-channel state. 
The upper figures are the predictions by fss2, while the lower ones by FSS. The ^Sq phase shift shows very strong 
Pauli repulsion, similar to the phase shift of the EN{I = 3/2) "^Si state. The EN{I ~ 1/2) "^Si phase shift predicted 
by fss2 starts from 180° at ps = 0, decreases moderately down to ~ 160° over the momentum region ps < 100 MeV/c. 
Then it suddenly decreases down to ~ 80° at around ps ~ 120 MeV/c. Beyond this momentum, a moderate decrease 
follows again. The situation in FSS is rather different. The ^Si phase shift predicted by FSS starts from 0° at ps = 0, 
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and shows a clear resonance behavior over the momentum range 100 < < 120 MeV/c. The peak value of the 
phase shift is about ~ 60° . In spite of the very different behavior of the diagonal phase shifts predicted by fss2 and 
FSS, the S'-matrices are found to be very similar to each other. Figure |l^ illustrates the Argand diagram showing the 
energy dependence of the 5'-matrix element ^33 — 77336^"^^^. Here i = 3 represents the 'EN{I — 1/2) ■^S'l channel. The 
resemblance of the circles given by fss2 and FSS is apparent. This implies that the strength of the central attraction 
in the EA^(/ — 1/2) channel is almost the same in fss2 and FSS. Figure |l^ also shows that the peak value of the 
RGM-H phase shift is about 45°, which indicates that the attraction of RGM-H is weaker than that of fss2 and FSS. 

The very strong reduction of 7^33 = |S'33| in Fig.|l^is related to the large enhancement of the transmission coefficients 
ryi3 and 7723 to the AA*" '^S'l {i — 1) and AA '^Di [i = 2) channels around ~ 120 MeV. These transmission coefficients 
are directly connected to the S~p An reaction cross sections (which we call process C), and the driving force for 
this transition is the one-pion tensor force. Table XI] shows contributions from each partial wave to this reaction cross 
section at ps = 160 MeV/c. The sum of the contributions from the transitions "^Si — > ^5*1 and ^Di amounts 

to about 120 mb both in fss2 and FSS, which is much larger than the contribution from the transition ^Sq 
(3.5 mb). This is in accordance with the analysis of the Ap system in the preceding subsection. Namely, there is a 
large cusp structure in the ■^S'l phase shift at the SA threshold, while a very small cusp is seen in the ^Sq channel. 
There is, however, some quantitative difference in the detailed feature of this tensor coupling between fss2 and FSS. 
In fss2, 7713 and 7723 are almost equal to each other over the momentum region where the experimental data exist 
(110 < < 160 MeV/c). This feature is very similar to RGM-F (see Fig. 6 in |^). On the other hand, 7713 < 7723 
holds in FSS, which is common with RGM-H (Fig. 5(b) in M). We can also read this difference from Table XII. The 
details of the cross section of the process C indicate that crpS'i — > '^Si) ^ a^^Si '^Di) in fss2, while a^^Si — > ■^S'l) 
< a{^Si ^Di) in FSS. 

For more detailed evaluation of S^p cross sections, it is important to take into account the pion-Coulomb correction. 
In [4^ , we incorporated the Coulomb force of the channel correctly in the particle basis, but the threshold energies 
of the £"71 and An channels were not treated properly. As is discussed in Sec. II D, we can now deal with the empirical 
threshold energies and the reduced masses in the RGM formalism, without spoiling the correct effect of the Pauli 
principle. Since the threshold energies and reduced masses are calculated from the baryon masses, these constitute a 
part of the pion-Coulomb correction in the YN interaction. Although the pion-Coulomb correction may not be the 
whole story of the CSB, it is certainly a first step to improve the accuracy of the model predictions calculated in 
the isospin basis. We can easily imagine that the small difference of threshold energies becomes important more and 



more for the low-energy S p scattering. In particular, the charge- exchange total cross section I] p 



(which 



we call process B) does not satisfy the correct law in the zero-energy limit, if the threshold energies of the S^p 
and E°7i channels are assumed to be equal. We therefore used the prescription [Q to multiply the factor {kf/ki), in 
order to get a(B) from (t(B) which is calculated by ignoring the difference of the threshold energies. Here ki and kf 
are the relative momentum in the initial and final states, respectively. We will see below that this prescription is not 
accurate and overestimate a(B) and cr(C). 

The largest effect of the pion-Coulomb correction appears in the calculation of the S^p inelastic capture ratio at 
rest rn [EtI- This observable is defined by [ESI 



1 



rR 



■ '"5=0 



with 



'^(5=0,1) (B) 
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where B and C denote the scattering processes S^p YPn and S^p — s- An, respectively. This quantity is the ratio of 
the production rates of the Yf^ and A particles when a particle is trapped in the atomic orbit of the hydrogen and 
interacts with the proton nucleus. For the accurate evaluation of rn, we first determine the effective range parameters 
of the low-energy S'-matrix using the multi-channel effective range theory. These are given in Tables XIII and XIV 
with respect to the ^6*0 and -I- ^Di states of the S~p - Y,^n - An system, respectively. The calculations are 
performed by using the particle basis with and without the Coulomb force. Figures |l^ and |l9| show the calculated 
phase shifts in the particle basis, when the Coulomb fo rce is incl uded. Also shown are the predictions by the effective 
range formula, using the parameters given in Tables XIII and XIV . For the ^Si + ^Di state, the effective range 
expansion breaks down around ps ~ 140 MeV/c due to the singularity of the matrix inversion. On the other hand, 
the effective range approximation works excellentl y for the ^Sq state. The scattering length matrices A {A'^) and A 
{A'^) are employed to calculate given in Table XV without (with) the Coulomb force. If we compare the result 
with the empirical values tr = 0.33 ± 0.05 |§, 0.474 ± 0.016 and 0.465 ± 0.011 H, we find that tr = 0.442 
predicted by fss2 is slightl y sm aller than the recent values between 0.45 and 0.49. The contribution from each spin 
state is also listed in Table XV. We find rs=o ~ 0.9, which indicates that the ac is very small in comparison with the 
CTB in the spin-singlet state. On the other hand, rs=i ^ 0.29 implies that most of the Y^N-AN channel coupling takes 
place in the spin-triplet state. We again find that the one-pion tensor force is very important in the T,N{I = 1/2) - 
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AN ^Si - ^Di channel-coupling problem. We also find that the effect of the Coulomb force plays a minor role for this 
ratio @. 

The inelastic capture ratio in flight rp, predicted by fss2 in the full calculation, is illustrated in Fig.|2C|. 

Unlike r^, this quantity is defined by using the total cross sections ctt(B) = (1/4)(To(B) + (3/4)cri(B) and (Tt(C) — 
(l/4)ao(C) + (3/4)ai(C) directly: 

aT{B) + aT{G) ■ ^ ' 

This is a rather sensitive quantity which depends on the relative magnitudes of the YPn and — > An total 

cross sections. In the momentum region of > 100 MeV/c, we find that rp in the particle basis gives rather similar 
values, irrespective of whether the Coulomb force is incorporated or not. The empirical value of r^? averaged over the 
momentum interval = HO - 160 MeV/c is rp = 0.47 ± 0.03 [|2|, and is plotted in Fig.|| by a cross. We find that 
the prediction of fss2, rp — 0.419, is too small, which is the same feature as observed in FSS [rp = 0.41) @. The 
main reason for this disagreement is that our S^p An cross sections are too large. 



F. Y N cross sections 



Figure displays the low-energy cross sections predicted by fss2 for the S~p and S+p scattering. The results 
of three different calculations are shown; the full calculation in the particle basis including the Coulomb force (solid 
curves), the calculation in the particle basis without the Coulomb force (dashed curves), and the calculation in the 
isospin basis (dotted curves). The effect of the correct threshold energies and the Coulomb force is summarized as 
follows. In the S+p scattering, the effect of the repulsive Coulomb force reduces the total cross sec tions by 11-6 
mb in the momentum range — 140 - 180 MeV/c, where the experimental data exist. (See Table XVI.) On the 
other hand, the attractive Coulomb force in the incident S~p channel increases all the cross sections. An important 
feature of the present calculation is the effect of correct threshold energy of the YPn channel. It certainly increases 
the S^p E'^n charge-exchange cross section, but the prescription to multiply the factor (kf/ki) overestimates this 
effect. The real increase is about 1/2 - 2/3 of this estimation. Furthermore, we find that this change is accompanied 
with the fairly large decrease of the elastic and S^p An reaction cross sections. Apparently, this effect is 
due to the conservation of the total flux. The net effect of the Coulomb and threshold energies becomes almost 
zero for the elastic scattering. The charge- exchange r eacti on cross section largely increases, and the T,~p — > An 
reaction cross section gains the moderate decrease. Table KVI summarizes this change of total cross sections in the 
momentum range = 110 - 200 MeV/c. We have also examined the effect of the threshold energies and the Coulomb 
force in a simple potential model, which flts the low-energy phase shifts of the model FSS. (The crosses shown in the 
FSS phase-shift curves in Figs.^(a), and 16 indicate the predictions of this potential model.) The cross section 
difference in this model is also shown in Table XVI for comparison. We flnd that both calculations give very similar 
results. Figure |2^ also shows the comparison with the experimental data. The flnal result given by fss2 reproduces 
the experimental data reasonably well, although the S^p An total reaction cross sections are somewhat too large. 

Figure G2^ shows the predicted differential cross sections by fss2 in the full calculation, compared with the experi- 



mental data. For E+p and elastic differential cross sections, the recent experimental data taken at KEK p4 55 
are also compared. The agreement between the calculation and the experiment is satisfactory. 

We show in Fig. |2^ the total cross sections of the NN and YN scatterings in the full energy range. The solid 
curves denote the fss2 result in the full calculation, while the dotted curves in the isospin basis. The "total" cross 
sections in the charged channels, pp, E^p, and S^p, are calculated by integrating the differential cross sections over 
the angles from cos6'min = 0.5 to cos6'max = —0.5. This is the reason the Coulomb result of the pp total cross sections 
at the higher energies Tiab > 400 MeV is very small, in comparison with the result in the isospin basis. Since the pp 
differential cross sections at these higher energies are very much V-shaped (see Fig. |^), the non-Coulomb calculation 
in the isospin basis is more reliable. The further difference between the dotted curve and the experiment in th pp 
scattering is due to the inelastic cross sections, which are zero in our single-channel calculation. The experimental 
analysis of the pp total inelastic cross sections shows that they are about 20 mb at Tiab = 800 MeV. For the np 
scattering, the inelastic contribution is smaller and is about 10 mb at the same energy. This implies that our fss2 
predicts the total elastic NN cross sections almost correctly up to the energies Tiab < 800 MeV. For the YN total 
cross sections, the pion-Coulomb correction is important only for the charged channels and the low-energy T,~p — > E^n 
reaction cross sections. 
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G. G-matrix calculation 



Figure shows saturation curves calculated for ordinary nuclear matter with the QTQ prescription as well as 
the continuous prescription for intermediate spectra. The results produced by the Paris potential ||5q| and the Bonn 
B potential are also shown for comparison. The ^-dependence of the nucleon, A and S s.p. potentials Usik) 
obtained with the continuous choice is shown in Fig.^ at three densities p — 0.5pQ, 0.7po and po, with pQ = 0.17 
fm""^ being the normal density. (These densities correspond to kp = 1.07, 1.2 and 1.35 fm~^, respectively.) For 
comparison, the results of the Nijmegen soft-core potential NSC89 |^8| calculated by Schulze et al. are also 
shown. The corresponding figures of the s.p. potentials predicted by our previous model FSS are given in Figs. 2 - 
5 of 1^. We find that fss2 gives a nucleon s.p. potential UN{k) very similar to that of FSS except for the higher 
momentum region qi > "i fm"^. As is discussed at the end of Sec. Ill A, the too attractive behavior of FSS in this 
momentum region is corrected in fss2, owing to the effect of the momentum-dependent Bryan-Scott terms involved 
in the S-meson and V-meson exchange EMEP. The saturation curve in Fig. |2j shows that this improvement of the 
s.p. potential in the high-momentum region has the favorable feature of moving the saturation density to the lower 
side, as long as the calculation is carried out in the continuous prescription. On the other hand, the saturation curve 
with the QTQ prescription suffers a rather large change in the transition from FSS to fss2. The prediction in fss2 
with the QTQ prescription is very similar to the prediction in Bonn model-B potential. It is interesting to note that 



our fss2 result is rather close to Bonn model-C for the deuteron properties (see Table [X), while to model-B for the 
nuclear saturation properties. The model-B has a weaker tensor force than model-C, which is a favorable feature for 
the nuclear saturation properties. 

We should keep in mind that the short-range part of our quark model is mainly described by the quark-exchange 
mechanism. The non-local character of this part is entirely different from the usual V-meson exchange picture in the 
standard meson-exchange models. In spite of this large difference the saturation point of our quark model does not 
deviate much from the Coester band, which indicates that our quark model has similar saturation properties with 
other realistic meson-exchange potentials. 

Figures p5|(b) and p5|(c) show the momentum dependence of the A and S s.p. potentials in nuclear matter obtained 
from the quark-model G-matrices of fss2. We find that the Uj^{qi) predicted by fss2 and FSS (shown in Fig. 3 of P" 
are again very similar for qi <2 fm~^. On the other hand, the repulsion of the Usiqi) predicted by FSS in Fig. 5 of"| 
is somewhat reduced. The partial wave contributions of the s.p. potentials Uf^jqi = 0) and U^{qi = 0) in symmetric 



nuclear matter at kp = 1.35 fm ^, predicted by fss2, are tabulated in Table XVII, together with the result of NSC89 



|39[| . The corresponding analysis in FSS is given in Table 1 of [|6[. For the A s.p. potential, the characteristic feature 
of fss2 appears in the less attractive ^Sq state and the more attractive '^Si state, in comparison with FSS. The partial 
wave contributions of fss2 now become very similar to those of NSC89, except for the "^Si + ^Di contribution. The 
extra attraction of fss2 to NSC89 in the A s.p. potential mainly comes from this channel (15 - 16 MeV). This is 
probably because the tensor coupling is stronger in fss2 than in NSC89. A minor excess of the attraction comes from 
the ^Pi + ^Pi and ^P2 + states (2 - 3 MeV). 

For the E s.p. potential, it should be noted that it is repulsive in the quark model, reflecting the characteristic 
repulsion in the '^Si + ^Di channel of the isospin / = 3/2 state (the Pauli repulsion). The repulsive feature of the 
S s.p. potential is supported by Dabrowski's analysis ||6^ of the recent {K^,Tr^) experimental data at BNL Jsi] . 
Quantitatively, the strength of the repulsion (which is 21 MeV in FSS for U^{qi = 0)) is reduced to 7 MeV in fss2. 
This change of the s.p. potential is mainly brought about by the 7 MeV reduction of the / = 3/2 -l-'^Di repulsion 
and the 4 MeV increase of the 1 = 1/2 "^Si + ^Di attraction. The latter feature is again related to the strong tensor 
coupling in fss2. On the other hand, the repulsive contribution of the / = 3/2 + ^Di state in NSC89 is very weak, 
since this channel has a broad resonance around = 500 - 800 MeV/c (see Fig. 1 in |42|). It is interesting to note 
that the attractive contributions to the A and E s.p. potentials from the 1 = 1/2 ^Si + '^Di state is more than 10 
MeV stronger in fss2, compared with those of NSC89. Although NSC89 is considered to be a model with a strong 
AN - 'SN coupling, the AiV - 'SN coupling of fss2 is even stronger. This feature should have some consequence in the 
energy spectra of the s-shell A-hypernuclei, if fss2 is used in the few-body calculations of these hypernuclei. 

The imaginary parts of the A and E s.p. potentials are shown in Fig.^(d) with respect to fss2. These results are 
rather similar to the predictions of FSS (see Fig. 4 in ref. |^). In particular, 3m U-s{qi = 0) at kp = 1.35 fm^^ is 
— 13.9 MeV in fss2 and —18.5 MeV in FSS. These results are in accord with the calculations by Schulze et al. ||5^ for 
NSC89. 

By using the G-matrix solution of fss2, we can calculate the Sheerbaum factor Sb, which represents strength of the 
s.p. spin-orbit potential defined through |10| 



Ui%r)=-^-Ss'-'-^i-a. (3.5) 
2 r ar 
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The explicit expression of SB{qi) (which actually contains the momentum dependence) in terms of the G-matrix is 
given in Eq. (50) of [|l^. Here we only consider Sb = Ssiqi = 0) as the measure of the s.p. spin-orbit strength 
in the bound states. The quark model description of the YN interaction contains the antisymmetric spin-orbit 
{LS^~^) component which originates from the FB LS interaction. The large cancellation between the LS and LS**^"^ 
contributions in the AA^ isospin 1—1/2 channel leads to a small s.p. spin-orbit potential for the A-hypernuclei. A 
very small ratio Sa/Sn < 1/10 was reported for FSS |jl^. In fss2 this cancellation is less prominent, since the present 
S-meson EMEP yields the ordinary LS component but no LS^~'> component (see Sec. II B). Since the total strength 
of the LS force is fixed in the A^A^ scattering, the FB contribution of the LS force is somewhat reduced. This can 
easily be seen from the simple formula given in Eq. (52) of [|lO|, which shows that in the Born approximation the FB 
LS contribution to the Scheerbaum factor is determined only by a single strength factor asx^niudC^b^ ■ The v alue of 
this factor is 29.35 MeV • fm^ for fss2, which is 3/5 of the value of FSS, 48.91 MeV • fm^. We show in Table pTvin 
the Scheerbaum factors Sb (at qi = 0) predicted by the G-matrix calculation of fss2 in the continuous prescription, 
with respect to the nuclear-matter densities p — 0.5po, 0.7po and pQ. At the normal density po with kp = 1.35 fm~^, 
we obtain for fss2 Sn = -42.4, Sa = -H.l and S^ = -23.3 (MeV • fm^), which gives the ratios Sa/Sn ~ 0.26 and 
Ss/Sn ~ 0.55. We find that Sa/Sn and S^/Sn become slightly smaller for lower densities. Each contribution from 
the LS and the iSl components in the even- and odd-parity states as well as / = 1/2 and / = 3/2 channels is 
shown in Table XIX for kp — 1.35 fm^^. The parenthesized numbers are the predictions by FSS. We find that a 
prominent difference between fss2 and FSS appears only in the contribution of Sa- Namely, the 5 MeV reduction 
of the LS"*^"-* contribution and the 2 MeV enhancement of the "^O LS contribution, which explains the increase 
of Sa = —3.5 in FSS to Sa = —11.1 in fss2. In the recent experiment at BNL, very small spin-orbit splitting is 
reported in the energy spectra of ^Be and |6^. A theoretical calculation of these AA^ spin-orbit splittings using 
OBEP AA'' interactions is carried out by Hiyama et al. The present result of fss2 is not entirely favorable for 
these experimental data. An accurate experimental determination of the s.p. spin-orbit strengths is very important 
to figure out the relative significance of the FB and EMEP contributions, both of which apparently constitute the 
bare two-body LS forces of the baryon-baryon interaction. 



IV. SUMMARY 



The purpose of this investigation is to construct a realistic model of the nucleon-nucleon (NN) and hyperon-nucleon 
(YN) interactions, which describes not only the baryon-baryon scattering quantitatively in the wide energy region, 
but also reproduces rich phenomena observed in few-baryon systems and various types of infinite nuclear matter. We 
believe that the present framework, incorporating both the quark and mesonic degrees of freedom into the model 
explicitly, is very versatile, since it is based on the natural picture that the quarks and gluons are the most economical 
ingredients in the short-range region, while the meson-exchange processes are dominant in the medium- and long- 
range part of the interaction. Since our quark model describes the short-range repulsion (which is observed in many 
channels of the baryon-baryon interactions) in terms of the non-locality of the quark exchange kernel, the effect of the 
short-range correlation is rather moderate, compared with the standard meson-exchange potentials. This can be seen 
in the magnitude of the Born amplitudes used in solving the Lippmann-Schwinger resonating-group equations (LS- 
RGM) [Ql and the Bethe-Goldstone equations ||^ , and also in the fairly reasonable reproduction of the single-particle 
(s.p.) spin-orbit strengths calculated in the Born approximation ||l0| . In [Q, we have seen that the Born amplitudes 
of the quark model have almost the same order of magnitude as the empirical scattering amplitudes obtained by 
solving the LS-RGM equation. The s.p. spin-orbit strength Sn predicted by the G-matrix solution of our quark 
model is almost equal to that in the Born approximation [ p^ , in contrast to the standard potential models like the 
Reid soft-core potential with the strong short-range repulsive core [|4j . Since the Born amplitudes in the quark model 
reflect rather faithfully the characteristic features of the LS-RGM solution, it is easy to find missing ingredients that 
impair the model. 

In this study we upgrade our previous model FSS ||^,^ in two respects. The first one is the renovation of the effective 
meson-exchange potentials (EMEP) acting between quarks. We extend our model to include not only the leading 
terms of the scalar and pseudo-scalar mesons but also the vector mesons with all possible standard terms usually 
used in the non-relativistic one-boson exchange potentials (OBEP). The second point is the exact incorporation of 
the pion-Coulomb correction in the particle basis. This includes the exact treatment of the threshold energies and the 
Coulomb exchange kernel, as well as the separate evaluation of the spin-flavor factors of the charged- and neutral-pion 
exchange EMEP. This improvement is necessary in order to study the effect of the charge symmetry breaking in 
the A^A^ and YN interactions. These two renovations require various mathematical techniques which are specifically 
developed in refs. and (l4j for these purposes. Appendix A in discusses a convenient transformation formula 
of the RGM kernel, which directly gives the Born kernel for the momentum-dependent EMEP at the quark level. A 



18 



procedure to avoid the problem of threshold energies in the RGM formalism is given in The new model fss2 

with these features has acquired much freedom to describe the NN and YN interactions more accurately than FSS. 
Three different types of calculations are carried out using fss2. The first one is the calculation in the isospin basis, 
which is used for determining the model parameters and also for the G-matrix calculation. The second and third 
calculations are performed in the particle basis with and without the Coulomb force. When the Coulomb force is 
included, the standard technique by Vincent and Phatak p9|] is employed to solve the Lippmann-Schwinger equation 
in the momentum representation. 

In the NN system, the incorporation of the momentum-dependent Bryan-Scott term and the vector-meson 
EMEP improves the quantitative agreement to the experimental data to a large extent. The momentum-dependent 
Bryan-Scott term, included in the scalar- and vector-meson EMEP, is favorable in extending our quark-model descrip- 
tion of the NN scattering at the non-relativistic energies to the higher energies up to Tiab = 800 MeV, and also in 
describing reasonable asymptotic behavior of the s.p. potentials in the high-momentum region. For vector mesons, we 
avoid the double-counting problem with the Fermi-Breit (FB) contribution by choosing small coupling constants 
around 1 especially for the flavor-singlet coupling constants ff'^ and fY"^- Since we have also chosen fg° = 0, the 
LS contribution from the vector mesons is very small. For the p- and i4r*-meson contributions, the selected value 
~ 2.6 through the parameter search is a standard size usually assumed in OBEP. Although the (//™)^ term 
usually gives the isovector spin-spin, tensor and quadratic spin-orbit (QLS) terms, we only retain the QLS term 
with the L^-type spin-spin term. This choice is rather ad hoc, but favorable since we do not want to introduce too 
strong cancellation between the one-pion tensor force and the p- meson tensor force in the '^S'l - '^Di coupling term 
of the NN interaction. Since the (3g) cluster wave function yields a large cut-off effect for the singular part of the 
one-pion exchange potential (OPEP), we also introduce a reduction factor cs for the spin-spin contact term of the 
OPEP central force, and multiply the short-range tensor term of the FB interaction by about factor 3. With these 
phenomenological ingredients, the accuracy of the model in the NN sector has now become almost comparable to that 
of the OBEP models. For the energies above the pion threshold, our single-channel calculation of the NN scattering 
seems to have given nearly satisfactory results, which are visible in the good reproduction of the differential cross 
sections up to Tiab = 800 MeV. The polarizations for the np and pp scattering have some unfavorable oscillations in 
the energy range Tjab = 400 - 800 MeV, but the improvement is a future work which definitely requires the explicit 
introduction of the inelastic channels such as the AiV channel. 

The existing low-energy data for the YN scattering is well reproduced. This includes; 1) Ap total cross sections 
at Pa < 300 MeV, 2) S+p and total and differential cross sections at pa < 200 MeV. The phase-shift difference 
of the ^Sq and '^Si states of the AA^ system at the maximum values is kept less than 10° in fss2, which seems to be 
necessary to describe the energy spectra of the hypertriton, and ^^^He systems. In the cusp region of the Ap total 
cross sections, the enhancement of the cross sections by the FB LS^~^ force is found in all versions of our quark model. 
The ^Pi resonance of the EA^(/ =1/2) state still remains in the original channel, which is a common feature found in 
both fss2 and RGM-H. The strong LS'^") force is one of the characteristics of the quark model, which is related to the 
spin-flavor SUe character of the AA^ - I]N{I = 1/2) channel coupling in the quark model. On the other hand, the ^S'l 
- '^Di coupling of the AA^ - UN {I — 1/2) system is caused by the one-pion tensor force. We find that the 5-matrix is 
very similar in fss2 and FSS, although the phase-shift behavior of the SAf(/ = 1/2) '^Si diagonal channel looks very 
different in these two models. We can conclude that the essential mechanism of the (5 + _D)-wave and P-wave AA^ - 
I]Af(/ = 1/2) channel couplings is unchanged among all versions of our quark model. 

We find a very small effect of the charge symmetry breaking for the Ap and An scattering, in contrast to the 
predictions by the Nijmegen group [Q. The energy region we are concerned with is pA < 300 MeV/c, which is too 
far away from the SA^ threshold to be affected by the small difference of the threshold energies ~ 3.5 MeV between 
E°n and T,~p, ~ 2.0 MeV between E+n and S^p, and the Coulomb attraction between S" and p. On the other 
hand, the low-energy T,~p scattering suffers a large effect of these threshold energies and the Coulomb attraction. 
In particular, we find that the effect of the correct threshold energies is very important for the detailed description 
of the low-energy T,~p total cross sections and the T,~p inelastic capture ratio at rest. The prescription to multiply 
(kf/ki) factor to reproduce the low-energy behavior of the Tt~p —> E^n charge-exchange total cross sections, used in 
our previous calculation in |^^, is not accurate enough to yield reliable estimates for the E~p elastic, — > E^n 
charge-exchange and E~p An reaction cross sections. In the final calculation in the particle basis with correct 
threshold energies, the increase of the E~p E^n reaction cross section is almost half of the (kf/ki) prescription in 
the momentum range < 200 MeV/c. Furthermore, the E~p elastic cross section and E~p — > An reaction cross 
section decrease fairly largely due to the effect of the flux conservation. The net effect of the pion-Coulomb correction 
on the T,~p elastic scattering is negligible. The E~p An reaction cross section seems to be still too large even in the 
present model fss2. This is reflected in the rather small values of the E~p inelastic capture ratios at rest {rn = 0.442) 
and in flight {rp = 0.419). 

The G-matrix calculation using fss2 shows that our previous results given by FSS are qualitatively pertinent. 
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In particular, the nucleon s.p. potentials in symmetric nuclear matter are very similar to the predictions of other 
realistic NN potentials. The nuclear saturation curve predicted by fss2 resembles the curve given by the Bonn model- 
B potential. It is interesting to note that the deuteron properties of fss2 are rather close to those of model-C, which is 
known to have a larger £>-state probability than model-B. Since fss2 reproduces the NN phase shifts at non-relativistic 
energies quite well, the difference of the off-shell effect between our quark model and the other OBEP models does 
not seem to appear so prominently, as far as the nuclear saturation curve is concerned. Some interesting features of 
our quark model appear in predictions for hyperon properties in nuclear medium. The A s.p. potential has a depth of 
about 48 MeV in the case of the continuous prescription for intermediate energy spectra, which is almost the same as 
the FSS prediction 46 MeV This value is slightly more attractive than the value expected from the experimental 
data of A-hypernuclei The S s.p. potential is repulsive, with the strength of about 7 MeV, which is smaller than 
21 MeV by FSS. The origin of this repulsion is the strong Pauli repulsion in the T^N{I = 3/2) '^Si state. This result 
seems to be consistent with the indication from the analysis by Dabrowski of the recent {K~ ,tt^) experiments |6^] 
at BNL. Future experiments will be expected to settle the problem of the E s.p. potential. One of the characteristic 
features of fss2 is the LS force generated from the scalar-meson EMEP. If this contribution is large, the cancellation 
of the LS and LS'^"-' components from the FB interaction becomes less prominent in the Scheerbaum factor Sa- The 
fss2 model predicts the relative ratio to N about Sa/Sn ~ 1/4, which is larger than the FSS value Sa/Sn ^ 1/12. 
The density dependence of the Sa/Sn ratio is rather weak in fss2. We should however keep in mind that this ratio 
is for the infinite N = Z system with the normal density po = 0.17 fm~^. We did not discuss in this paper the most 
appropriate relative strength of the LS terms which come from the FB interaction and the scalar-meson EMEP, since 
this cannot be determined only from the NN data. We definitely need more experimental information concerning 
each contribution of the LS and LS'^~^ forces in the YN interaction. 

Finally we note that it is an important future subject to consider few-body systems including hyperons in the 
scope of the quark-model baryon-baryon interactions. The hypertriton calculation can be performed in the Faddeev 
formalism and the stochastic variational approach by using the quark-exchange kernel directly for the NN and YN 
interactions. The study of hyperonic nuclear matter is also interesting, since the G-matrix calculation of AA and 
SN interactions ||6^ in the models FSS and fss2 is now in progress. Since the E s.p. potential is repulsive in the 
quark-model description, the admixture of the E particle is suppressed, and this should affect the behavior of the A 
particles in dense hyperonic nuclear matter. 



ACKNOWLEDGMENTS 



This research is supported by Japan Grant-in- Aid for Scientific Research from the Ministry of Education, Science, 
Sports and Culture (12640265). 



APPENDIX A: EMEP EXCHANGE KERNEL 

In this appendix we extend the derivation of the EMEP exchange kernel developed in Appendices A and B in [Q, 
to deal with various interaction pieces of the V mesons, including the LS and QLS terms. The Coulomb exchange 
kernel and internal-energy contribution from EMEP are also discussed. 

The systematic evaluation of the quark-exchange kernel is carried out by assuming a two-body interaction 

(7.,=^aX'<', (Al) 
n 

where represents the spin- flavor part (the color part is wfj — 1 for EMEP) and ufj the spatial part. Four 
different types of the spin-flavor factors fl = C, SS, T, LS are required for the most general EMEP up to the V 
mesons; w'-^ = 1, w^^ = (cri • 0-2), = [(7i(72]|f\ and w'"^ = {cri -I- (T2)/2. For the flavor octet mesons, these spin 
operators should be multiplied with (A^Aj), where A represents the Gell-Mann matrix in the flavor SU^ space. The 
spin-flavor factors are defined by Eq. (A. 3) of J^lfor each w^j with the quark exchange number x = 0, 1 and 
the five interaction types T = E, S, S", Z?+, D_ |67|. The non-central factors are defined by the reduced matrix 
elements for the tensor operators of rank 1 and 2. For example, the tensor operator is expressed as 

2 



Si2{k, fe) = 3 (cTi • k) {(T2 ■ k) - (cTi • (T2) k 

, — r (o\ 1 (0) 

= 3\/l0 [cjixa2r^y2{k) , (A2) 
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where y2^,{k) = yi^TSfe^Fa^ (fc 



The reduced matrix elements of the spin operators at the baryon level are 

assu med to be 1. For the spatial part, we also need three extra types = C(l), SS{\), QLS listed in Table 
XX. This table shows the polynomial functions u{k,q) accompanied with the Yukawa function in the momentum 
representation through 



u{k,q) = 



An 



■u{k,q) , 



(A3) 



and the spatial part of the Born kernel M[^(<7/, defined in Eq. (A. 4) of explicitly. The formulae Eqs. (A. 18) - 
(A. 21) given in greatly simplify the procedure to obtain these results. The spatial functions f^{6) are explicitly 
given below. 

In Eq. (Al) the coefficients and the correspondence among il, fl' and fl" are tabulated in Table XXl. The 



(A4) 



EMEP contribution of the exchange Born kernel in Eq. ( |2.4| ) is calculated through 

M'\qf,q.) 0'Hqf,q,) = J] X^^ M^^ {qj, q,) . 

T 

The final result is as follows. For the central part, we have — C, C(l), SS, SS{1) types with 



M^(v)(g^,q,) 



-9 



T 
■,2 



M^^'r.)i^q,,q,)^2,-{\Wx?rfr''\0), 



M^^(v)(q^,q,) 



f2l I rn 



f2 2 

-fib' 



T 



[0) 



(A5) 



Here 7 — {m/2mud) and cos9 = {qf ■ In these central terms, the spin-flavor factors xf^^^ should be replaced 

C SS 

with —X-^^g, , because of the subtraction of the internal-energy contribution in the prior form. The tensor parts of 
the PS and V mesons are given by 

M^("-)iqf,q,) = 



(A6) 



where the V- meson contribution is also given for completeness although this term is not used in fss2. The EMEP 
QLS contribution reads 



M'?^^(v)(g,,5,)=(^ 

but also contains the tensor contribution 

MQ^i^){qf,q, 



-fib' 



X 



^n,fS','{0)~X^^_fS'Ji9) 



9'b' 

_f2 8 2 

J m 3 / 



Am? 



(A7) 



(A8) 



which we call — QT term. In Eqs. (A7) and (A8), the QLS contribution from the S meson is also shown for 
completeness, although this term is negligibly small in fss2. The LS term has the contribution both from the S meson 
and the V meson: 



M^'5(v)(q,^.,q,,) 



(A9) 
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For the tensor and QLS tensor terms in Eqs. (A6) and (A8), each interaction term with T = S, S", D+, D- types 
should be rearranged to = T, T' , T" types in Eq. (2^), according to the rules given in Eq. (B.13) or (B.17) of [Q. 

The EMEP spatial functions f^{e) used here are defined by extending /^^(6I), and /J^(6l) given in 

Eq. (B.18) of 1^]. The following four basic functions are used in Table XX. 



exp{-i62(q2 + fc2)| 55^^(0) 

(i)*exp{-i62(q2 + ife2)}5^.,^ {\h\k\) 
(I)*exp{~i62fc2}5^.._ (^%l) 



for r = <^ 



E 
S 

D+ 



fT{o) = f^{e) with 3^„,(p) 



/l^(e) = /r(^) with y^Ap) 



2aT 



-47r 



: A)i exp {- [4(g2 + fc^) _ fe . g] } zf^ = + fe| 



exp{-i62(g2 + ifc2)} (i6|fe|) 



for r 



(D^xpl-ifc^fc^} Zi'^^ (^6|g| 



r ^ 



(AlO) 



The 5"-type spatial function f^g>{0) is obtained from fg{0) by taking k k. There is no £^-type possible for the 
non-central terms. The coefficients ar are given by as — as' = (ll/8)a£;, ao^ — '2aE, and a£)_ — {3/2)aE, with 
as = (m6)^/2 = (l/2)(mc6/?i)^. For the spin-spin part of the one-pion exchange EMEP, ya-rip) — (l/2ar) should be 
modified into yarip) ~ (l/2ar)- The modified Yukawa functions ya{p), 2a \p), and Za{p) are essentially given 
by the error function of the imaginary argument: 



ya{p) = 

2i'\p) 



2a 



dt . 



(All) 



The other spatial functions appearing in Eqs. ( A5) - ( |A9| ) are defined by using the four spatial functions in Eq. ( A1C| ): 



8aE 



> + 



1 

2m 



Q 

fe2 





jik + q) 





} for r = <; 



r E 

s 

D+ 



''Note that f^{e) = ~fr^{e) and /f°(6') = ^ fr^ {d) except for the difference of but /r"(6l) here contains different 
numerical factors from those of fr^ (9). 
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2m2 



+7r^ < 



-(- 

4 V m 



V for r = <; 



r 

5 



/r(^) = <! (^) - for r = <; 



(A12) 



For the numerical calc ula tion s, it is convenient to include the direct term also in the above expressions. This can 
be achieved in Eqs. (^) - (^), if we further add X^^^^f^" (0) term, in addition to the X%^fj^'^{9) term. The 

direct- type spatial functions f^{0) are given by 



47r 



fe2 



re 3^ 



(A13) 



and 



-1 

1 

3aE 

-(-V 

9 Vm/ 

]_ 

3aE 



for = < 



CD 

C(l) 
55(1) 

gL5 



(A14) 



The Coulomb exchange kernel is very similar to the color-Coulombic term of the FB interaction, as is discussed in 
Sec. II D. Only difference is 1) as ^ a = {e'^/fic), 2) the definition of the Coulomb spin-flavor factor 



T 



(A15) 



and 3) the modification of the spatial function ho{p) in Eq. (B.5) of 0, by the effect of the Coulomb cut-off at Re- 
The last modification is achieved by 



ho{p) ho{p) - g{x,p) 
g{x,p) = e" 



cos{2 pxt) dt 



(A16) 







with X = (l/V2)(i?c/6), (2/VTT)(i?c/&), {ll2){Rclh) and {1/ ^/i){Rc /h) for the T ^ E, S or 5', D+ and D. 
types, respectively. The function g{x,p) is expressed as 
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g{x,p) = e 5to erf (a; + ?p) 

2p 

= y^e~'^ [sm{2px) '^e w{p + ix) + cos{2px) '<sm w{p + ix)] , (A-17) 
2p 

where w{z) — e^^ erfc {—iz) with erfc (2;) = 1 — erf (z). We note the simple relationship 

5(0. P) = 1^ '^m w{p) = ho{p) , 
2p 

g{x,Q) = e--" . (A18) 

For example, the T = E type spatial function is given by f^^{0) ~ ^2/7rQ!a;TO„rfC^(4/3)/(0) ^1 — e"'^^/*"^^/^^ , since 
^(0) = 1 (Cf. Eq. (III)). 



The EMEP contribution to the internal energies of the octet baryons originates only from the central force. It reads 
El, = mg^ ( -1 + ^) r„,(0) - mg^ X^i ^ Y^JO) , 

En = ^f(^Xxo'EYZ{0), 
3 \m^+ J 

El, = mf,' X^^ {l + ^ (0) + mfj X^ ^ (l + C (0) , (A19) 

where the values of the modified Yukawa functions at the origin are given by Ya{Q) = l/y^Tra — e"efrc(Y^) and 



(0) = Yq(0) — l/(2Q!0ra). For the qq interaction in Eq. (2.11), the last spin-spin terms of Ef^, and iJj^t should 
be modified into and mfm^ X^^ {2^'^/iaE) ^^^(0), respectively. 

APPENDIX B: SPIN-FLAVOR FACTORS IN THE EMEP SECTOR 

All the spin- flavor-color factors in the quark sector are already published in refs. [Q, 0, and jl^. Some simple 
spin-flavor factors including the direct factors in the EMEP sector are discussed in our previous publication [H, 
together with the details of the singlet-octet meson mixing]^ Since the V mesons and the non-central terms of the S 
mesons are not introduced in the LS factors are not shown. Here we show these LS factors and a complete list 
of the spin- flavor factors in the isospin basis with respect to the exchange terms in the EMEP sector. The Coulomb 
exchange factors are also shown only for the NN system. The detailed derivation of these factors and the spin-flavor 
factors in the particle basis (tt* - 7r° separation), as well as the exchange Coulomb factors for the YN system, will 
be pubhshed elsewhere. 

The LS factors for the flavor singlet and octet mesons are expressed by the electric-type {^Ixx) ) ^'^'^ the magnetic- 
type (e^xX)) vectors as follows (we use the same notation as in Appendix of [Q): 

[ (AA) ^ (00) ] 

K 
X[ 



L5(00) 
OB 




LS(OO) 
IE 


= Xn{S=1)- 




L5(00) 
15 


^Xn{S^\) + 




L5(00) 
1D+ 


— ^^IB > 




LS(OO) 
ID- 


^-^(e'^tpm^ 

12 ^ 


e^te^) 



'A different notation = CN, TN is used in Appendix of M|, to specify Q — C, T used here. 
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where Xn{S = 1) = -(1/12) (e'=te'= + e'"1'e™). We have simple relationships Xq^'-°°'^ +^0^0^°^ = 5 and 
(AA) = (11) 



Er^fr'""^ = 5XAr(^ = l) 



T^LS(ll) 


-v^-LS(ll) 
~ 3 ' ^0D+ 






4 „LS(00) 

- 3^ii. , 






4 „LS(00) 
- 3^15 . 








I) 




-2 





(B2) 

The flavor octet factors X^^^^ are further separated into the sum of the isoscalar, isovector, and 1—1/2 (strange) 
factors in the isospin representation. This separation is achieved in two steps. First we calculate the reduced matrix 
elements 

X?rW = {-9)UPs,^'^\\j2wf/^U'^) , (B3) 
3 ^ — / J miTOo 

with respect to wff = 1, (cr, • a-j) , (cr, +crj)/2 (1||5'||1)"\ and [a^ x CTj]^^^ (1|| [as^ x ctsJ^^^ ||1)"^ Here the SU3 
symmetry breaking is introduced by the parameter A through 



mud 



Y, . (B4) 



The flavor-singlet factors are obtained by setting A = 1: X^q-i^) = ^^7"°^- the other hand, X^^(— 1/2) is the 
contribution from the As vertex. For the product of two coupling constant operators (/i cos 9 + fs sin OX^i) (/i cos 9 + 
/s sin 9Xsj ) , the isoscalar spin- flavor factors including the (quark) coupling constants are obtained through 

W^ri^) = i (\/3/i cos 9 + fs sin 0) ^,"^(1) Vs/i cos 9 + X^^ ^ fs sin , (B5) 

for the T ^ E types. (We do not need X^^, since it is replaced by —X^g, due to the subtraction of the internal 
energy part.) The other isoscalar factors are obtained from W^'j-{9 — tt/2). The non-isoscalar part of the spin-flavor 
factors are obtained from 

X^r'"=xT'-lx2r{-l) . (B6) 

(Note the hypercharge operator is given by y = As /Vs.) These are further separated into the isovector and 7=1/2 
factors by using the standard recoupling techniques of the SU3 algebra. We will specify this isospin dependence by 
using names of the S mesons 6 and k for = C, LS, and those of PS mesons tt and K lor ^1 = SS, T. These factors, 
however, can also be used for the vector mesons. 



In the following, X[^(A) in Eq. (B3) is expressed as ^[7-- The spin and isospin operators, cr^ and r^, are with 



respect to the two baryons Bi with i = 1 or 2 in the initial state, and the flavor exchange operator Pp is supposed to 
operate on the ket state. For example, (xi • T2) = /(/-f 1) — 11/4 for the EA^ channel. We assume B2 = B4 = N, unless 
otherwise specified. The spin-fiavor-color factor for the exchange normalization kernel is given by Xn = X^j^^ (A — 1). 
The factors of the interaction type T — S' are all equal to the T — S type factors for the diagonal channels, since they 
are Hermitian conjugate to each other. For the KN - TjN coupling terms, the S type and the S' type are interchanged 
if the initial state and the final state are interchanged. (See § 5 of for the symmetries of the spin-fiavor factors.) 



For f2 = LS, these factors with T = S, S' are not shown since they are not necessary, as in Eq. (A9) 
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1. Spin-Flavor factors of Eq. (B3) 
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(B7) 
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APPENDIX C: DEUTERON WAVE FUNCTIONS 

The relative wave functions for the deuteron in the momentum representation, /^(q) ^ 1/(7^ + <f')Ti-{q, ~*7, ^^d), 
satisfy the homogeneous equation 



(CI) 



where Vu'{p, q, —td) is the partial-wave components of Eq. (p^). Since /^(g) are the relative wave functions of the 
RGM equation, one needs to renormalize them through the square root of the normalization kernel This can be 
achieved by calculating 



'X^^ and X^]^ factors in Eq. (A2) of M] contain misprints in the (cti ■ (T2) terms. 
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FeXq)=qffXq) + qY,RNe {q,b^/3) 



In 



N 



V^ + lN + 1 



'Ni 



(C2) 



where RNeif, v) represents the radial part of the harmonic-oscillator wave function with the width parameter and 
In = (1/3)^"*"^ with iV = 0, 2, 4, • • • are the eigen- values of the exchange normalization kernel for the states of 
the NN system. The harmonic-oscillator components J^v^ of fi{q) are calculated from 



JNt= I q^dq RNt{q,b^ /i) fi{q) 



(C3) 



The deuteron wave functions ui(r) in the coordinate representation (customarily written as u{r) = uo{r) and w{r) = 
U2(r) for the 5*- wave and D-wave states, respectively) are obtained from the Fourier transformation 



ue{r) = / 

In particular, fe{q) are normalized such that 



dq {.qr)jg{qr) Ft{q) 



Jo p Jo 



{1)f 



(C4) 



(C5) 



We follow the standard ansatz |69 3l|,|ll| for the simple parameterization of the deuteron wave functions: 



TT q- + 



for 



£^0 

i = 2 



-7jr 



ue{r) 



for 



e = o 
e^2 



(C6) 



The range parameters 7^ are chosen as 7^ = 7 + (j — 1)70 with 70 = 0.9 fm ^ and n = 11. The coefficients Cj {j = 1 
- 10) and Dj {j = 1 - 8) with 7 = 0.23186542 fm~^ are given in Table XXII for the deuteron wave functions in the 
full calculation. The other coefficients, namely, the last Cj and the last three Dj should be calculated from Eqs. (C.7) 
and (C.8) of M. 
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TABLE I. Quark-model parameters, SU3 parameters of the EMEP, S-meson masses, and some reduction factors cg etc. 
for the models fss2 and ESS. The p meson in fss2 is treated in the two-pole approximation, for which mi and m2 (/32) are 
shown below the table. 





b (fm) 


m„d (MeV/c^) 


as 


A = ms/mud 


fss2 


0.5562 


400 


1.9759 


1.5512 


ESS 


0.616 


360 


2.1742 


1.526 




f! 


Is 




e| ^) 


fss2 


3.48002 


0.94459 


33.3295° 


55.826° 


ESS 


2.89138 


1.07509 


27.78° 


65° 






rPS 
78 






fss2 


- 


0.26748 


- 


(no r;,??') 


ESS 


0.21426 


0.26994 


-23° 






rVc 
.11 


78 


rVm 
7l 


rVm 2) 
78 


fss2 


1.050 





1.000 


2.577 


(MeV/c ) 


rric 


mg. 


ms 




fss2 


800 


1250 


846 ^> 


936 


ESS 


800 


1250 


970 


1145 




cs 


Cqss 






fss2 


0.4756''' 


0.6352 


3.139 




ESS 


0.381 









1) el is used only for S7V(/ = 3/2). 

2) ff^ — 35.264° (ideal mixing) and two-pole p meson with mi (/3i) 
912.772 MeV/c^ (0.48747) (2| are used. 

3) For the NN system, = 720 MeV/c^ is used. 

4) Only for tt, otherwise 1. 

5) The enhancement factor for the Fermi-Breit tensor term. 



664.56 MeV/c2 (0.34687) and ma [P2) 



TABLE IL The interaction types and the meson species introduced in the EMEP of the models fss2 and ESS. C represents 
the central force, SS the spin-spin force, T the tensor force, and QLS the quadratic spin-orbit force. C{BS) implies that the 
momentum-dependent Bryan-Scott term is also included for the central force. The tensor term of the vector mesons is switched 
off. 



model 


meson type 


interaction type 


mesons 




S 


C{BS) + LS 


e, S*, 5, K 


fss2 


PS 


SS + T 


TT, K 




V 


C{BS) + LS + QLS 


uj, (j), p, K* 


ESS 


s 


C 


€, S* , S, K, 




PS 


SS + T 


77', r], TV, K 
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TABLE III. Quark and EMEP contributions to the N - A mass difference (A_Ejv-a) and the A - E mass difference ( A_Ea-e ) 
in MeV, calculated in the isospin basis. The model is fss2. The mass ratio of strange to up-down quarks, A = {ms/rriud) = 1.5512, 
is employed to calculate the quark contribution in A_Ea-e . The details of the EMEP contribution to A_Ea-e are given in Table 
. See Table [fl for the two-pole p-meson parameters and the other EMEP parameters. 





P 


(MeV/c^) 


E (McV) 




Quark 




293.33 




S 


720 


-164.70 


AEn-a 




138.039 


71.56 




OJ 


781.940 


-34.36 






1019.413 


-0.19 




P 


two-pole 


80.59 


exp 


293.3 


total 


246.23 


ASa-s 


Quark 


(A = 1.5512) 


69.49 








7 QS 


GXp 


77 AA 


lOltii 


77 /1 7 




TABLE IV. The baryon 


masses employed in the particle-basis calculation l3(T| . 

1 tj tr ir__jj 




B 






Mb (MeV/c^) 


P 






938.2723 


n 






939.565 


A 






1115.683 


E+ 






1189.37 


EO 






1192.642 


E" 






1197.449 
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TABLE V. Details of the EMEP contributions to the baryon mass difference (in MeV) in the isospin and particle bases. 
In the particle basis, only the mass difference of the charged and neutral pions is introduced. 



13 


(MeV/c^) 


isospin basis 


particle basis 




ASa-e 


A-Ea_eo 




5 


846 


-87.345 


-87.345 





K 


936 


75.072 


75.072 





■K 


138.039 


47.704 


- 


- 




139.570 


- 


23.747 


-11.873 




134.976 


- 


24.061 


12.031 


K 


495.675 


-32.716 


-32.716 





U) 


781.940 


-26.587 


-26.587 





4> 


1019.413 


-1.799 


-1.799 





p 


two-pole 


53.723 


53.723 





K* 


893.880 


-20.072 


-20.072 





total 




7.980 


8.084 


0.158 



TABLE VI. Calculated threshold energies (in MeV) compared with the empirical values in the non-relativistic kinematics. 
Note that the effect of the charged and neutral pion mass difference on the n and p internal energies is zero. 



channel 


eQ 




ECi^ 




riexp 
th 


AiV 
EiV 




69.488 




7.980 








77.468 




77.471 


Ap 

E+n 




69.488 
69.488 




7.926 
8.084 









77.415 
77.572 




74.980 
76.959 


An 
E°n 
E-p 




69.488 
69.488 




8.084 
7.926 





2.066 




77.572 
79.480 




76.959 
80.473 


E-p 
E°n 
An 





-69.488 




0.158 
-7.926 




-2.066 
-2.066 




-1.908 
-79.480 




-3.514 
-80.473 
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TABLE VII. The threshold energies in the relativistic kinematics, used in the particle-basis calculation. The momentum 
is measured in MeV/c and the energy in MeV. The difference, AS'°* = iJ^^^ (relativistic) — Sj^', is given in the last column. 
Here B^^^ (relativistic) = (pt™)^/2/Ltinc with /Uinc being the non-relativistic reduced mass of the incident channel. 



channel 




relativistic 




non-relativistic 




dif 




Pth 


Pth 






iTical 
Eth 


AS'"* 


Ap 

















— 


S+n 


279.04 


633.14 


76.389 


74.980 


77.415 


-1.026 


S°p 


282.77 


642.20 


78.443 


76.959 


77.572 


0.871 


An 



















E°n 


282.87 


641.93 


78.441 


76.959 


77.572 


0.869 


E-p 


289.38 


657.79 


82.093 


80.473 


79.480 


2.613 


E-p 



















E°n 






-3.652 


-3.514 


-1.908 


-1.744 


An 






-82.093 


-80.473 


-79.480 


-2.613 
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TABLE VIII. The np phase-shift parameters calculated in the isospin basis (in degrees). The results given by OBEP, Paris 
and Bonn potentials are cited from Table 5.2 in ||3l|. 



State 


Model 






Tlab 


(MeV) 










25 


50 


100 


150 


200 


300 




fss2 


80.98 


63.03 


43.21 


30.51 


21.00 


7.02 




OBEP 


80.32 


62.16 


41.99 


28.96 


19.04 


4.07 




Paris 


80.35 


62.28 


42.26 


29.24 


19.25 


3.91 




Bonn 


80.30 


62.19 


42.27 


29.64 


20.31 


7.06 




SP99 


80.26 


62.10 


42.22 


29.69 


20.51 


7.07 




fss2 


-2.82 


-6.52 


-12.43 


-16.59 


-19.49 


-22.58 




OBEP 


-2.99 


-6.86 


-12.98 


-17.28 


-20.28 


-23.72 




Paris 


-2.95 


-6.77 


-12.85 


-17.22 


-20.42 


-24.52 




Bonn 


-3.03 


-6.98 


-13.25 


-17.64 


-20.62 


-23.43 




SP99 


-2.72 


-6.84 


-13.09 


-16.69 


-19.08 


-23.04 




fss2 


1.68 


1.91 


2.21 


2.68 


3.33 


4.97 




OBEP 


1.76 


2.00 


2.24 


2.58 


3.03 


4.03 


ei 


Paris 


1.69 


1.89 


2.14 


2.59 


3.21 


4.76 




Bonn 


1.82 


2.08 


2.29 


2.54 


2.82 


3.19 




SP99 


1.69 


2.14 


2.91 


3.55 


4.08 


5.06 




fss2 


-6.70 


-10.26 


-14.82 


-18.38 


-21.57 


-27.32 




OBEP 


-7.21 


-11.15 


-16.31 


-20.21 


-23.47 


-28.70 


'Pi 


Paris 


-7.11 


-10.95 


-15.72 


-19.08 


-21.73 


-25.92 




Bonn 


-6.90 


-10.48 


-15.11 


-18.88 


-22.41 


-29.17 




SP99 


-6.71 


-9.98 


-14.47 


-18.29 


-21.56 


-26.62 




fss2 


3.67 


8.82 


17.09 


22.26 


25.06 


26.38 




OBEP 


3.88 


9.29 


17.67 


22.57 


24.94 


25.36 




Paris 


3.96 


9.60 


18.64 


24.19 


27.15 


28.54 




Bonn 


3.88 


9.27 


17.41 


21.68 


23.09 


20.84 




SP99 


3.87 


9.37 


17.89 


22.73 


25.03 


25.47 




fss2 


52.26 


41.94 


27.51 


16.91 


8.41 


-4.86 




OBEP 


50.72 


39.98 


25.19 


14.38 


5.66 


-8.18 


'So 


Paris 


48.38 


38.12 


23.88 


13.38 


4.84 


-9.01 




Bonn 


50.03 


39.15 


24.36 


13.72 


5.30 


-7.62 




SP99 


51.30 


41.88 


28.24 


16.95 


7.74 


-5.49 




fss2 


8.55 


11.25 


9.04 


4.02 


-1.49 


-12.10 




OBEP 


9.34 


12.24 


9.80 


4.57 


-1.02 


-11.48 




Paris 


9.21 


11.93 


9.83 


5.32 


0.48 


-8.52 




Bonn 


9.57 


12.79 


10.88 


6.02 


0.66 


-9.66 




SP99 


8.24 


10.75 


8.18 


3.15 


-1.95 


-11.63 




fss2 


-5.23 


-8.68 


-13.45 


-17.27 


-20.77 


-27.26 




OBEP 


-5.33 


-8.77 


-13.47 


-17.18 


-20.49 


-26.38 


'Pi 


Paris 


-5.27 


-8.64 


-13.44 


-17.35 


-20.91 


-27.30 




Bonn 


-5.17 


-8.53 


-13.38 


-17.62 


-21.73 


-29.87 




SP99 


-4.75 


-8.15 


-13.52 


-17.92 


-21.64 


-28.06 




fss2 


0.64 


1.47 


3.29 


5.30 


7.27 


10.28 




OBEP 


0.68 


1.58 


3.34 


4.94 


6.21 


7.49 


'D, 


Paris 


0.78 


1.85 


4.00 


5.90 


7.47 


9.19 




Bonn 


0.72 


1.72 


3.76 


5.62 


7.04 


8.32 




SP99 


0.64 


1.59 


3.60 


5.60 


7.33 


9.75 
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TABLE VIII. -continued 



State 


Model 






Tiab (MeV) 












zo 


OU 


1 nn 
100 


ioU 


zUO 


300 




fss2 


2.58 


6.26 


12.43 


15.92 


17.36 


16.97 




OBEP 


2.62 


6.14 


11.73 


14.99 


16.65 


17.40 


3 D 

-P2 


Paris 


2.61 


5.97 


11.34 


14.68 


16.39 


16.74 




Bonn 


2.54 


5.89 


11.14 


14.24 


15.93 


17.22 




bFyy 


2.70 


5.93 


10.92 


14.11 


16.05 


17.83 




fss2 


0.10 


0.32 


0.72 


0.98 


1.08 


0.75 




OBEP 


0.11 


0.34 


0.77 


1.04 


1.10 


0.52 


3 ITT 

F2 


Paris 


0.11 


0.36 


0.79 


1.05 


1.06 


0.49 




Bonn 


0.11 


0.35 


0.81 


1.14 


1.28 


0.87 




SP99 


0.09 


0.33 


0.85 


1.19 


1.31 


0.90 




fss2 


-0.82 


-1.77 


-2.85 


-3.22 


-3.24 


-2.94 




OBEP 


-0.86 


-1.82 


-2.84 


-3.05 


-2.85 


-2.02 


£2 


Paris 


-0.87 


-1.80 


-2.73 


-2.90 


-2.74 


-2.14 




Bonn 


-0.85 


-1.77 


-2.74 


-2.97 


-2.84 


-2.23 




SP99 


-0.70 


-1.48 


-2.37 


-2.71 


-2.74 


-2.29 



TABLE IX. Deuteron properties by fss2 in three different calculational schemes, compared with th predictions of the Bonn 
model-C potential [Bl| and the experiment. 



isospin basis particle basis Bonn C Expt. Ref. 







Coulomb off 


Coulomb on 






td (MeV) 


2.2250 


2.2261 


2.2309 


fitted 


2.224644 ± 0.000046 ||] 


Pd (%) 


5.490 


5.490 


5.494 


5.60 




77 = Ad /As 


0.02527 


0.02527 


0.02531 


0.0266 


0.0256 ± 0.0004 || 


rms (fm) 


1.9598 


1.9599 


1.9582 


1.968 


1.9635 ± 0.0046 ||] 


Qd (fm^) 


0.2696 


0.2696 


0.2694 


0.2814 


0.2860 ± 0.0015 Hi 




0.8485 


0.8485 


0.8485 


0.8479 


0.85742 
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TABLE X. Effective range parameters of fss2 for the A''A^ and YN interactions in the single-channel formula. For the pp and 
nn systems, the calculation in the particle basis uses /f x 0.9949, in order to incorporate the effect of the charge independence 
breaking. Unit of length is in fm^^+^ in a, fm"^^"*"^ in r and fm~^^ in P for the partial wave £. The experimental values are 
taken from ||, @, g], Q, for iViV, ^ for E+p, and H, H for Ap. 







isospin basis 


particle basis 
Coulomb off 


Coulomb on 


Expt. 


pp 'So 


a 
r 
P 


-23.76 
2.584 
0.0393 


-17.80 
2.675 
0.0423 


-7.810 
2.574 
0.0334 


-7.8063 ± 0.0026 
2.794 ± 0.0014 


PP ^Po 


a 
r 
P 


-2.740 
3.867 
-0.014 


-2.876 
3.831 
-0.0130 


-3.004 
3.312 
-0.0125 


-4.82 ± 1.11, -2.71 ±0.34 
7.14 ±0.93, 3.8 ± 1.1 


pp ^Pl 


a 
r 
P 


1.740 
-8.196 
0.0009 


1.821 
-8.159 
0.0010 


2.112 
-8.269 
-0.0063 


1.78 ±0.10, 1.97 ±0.09 
-7.85 ± 0.52, -8.27 ±0.37 


nn 'So 


a 
r 
P 


-23.76 
2.584 
0.0393 


-18.04 
2.672 
0.0423 


-18.05 
2.672 
0.0423 


-18.5 ±0.3, -18.9 ±0.4 
2.75 ±0.11 


nn ^Po 


a 
r 
P 


-2.740 
3.867 
-0.0140 


-2.881 
3.823 
-0.0131 


-2.881 
3.822 
-0.0131 




nn ^Pi 


a 
r 
P 


1.740 
-8.196 
0.0009 


1.823 
-8.151 
0.0010 


1.823 
-8.152 
0.0010 




np 'So 


a 
r 
P 


-23.76 
2.584 
0.0393 


-27.38 
2.528 
0.0324 


-27.87 
2.525 
0.0324 


-23.748 ±0.010 
2.75 ± 0.05 


np ^Po 


a 
r 
P 


-2.740 
3.867 
-0.0140 


-2.466 
3.929 
-0.0186 


-2.466 
3.929 
-0.0186 




np ^Si 


a 
r 
P 


5.399 
1.730 
-0.010 


5.400 
1.730 
-0.0096 


5.395 
1.730 
-0.0097 


5.424 ± 0.004 
1.759 ± 0.005 


np 'Pi 


a 
r 
P 


2.824 
-6.294 
-0.0058 


2.826 
-6.299 
-0.0058 


2.826 
-6.299 
-0.0058 




np ^Pi 


a 
r 
P 


1.740 
-8.196 
0.0009 


1.582 
-8.185 
0.0004 


1.582 
-8.185 
0.0004 
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TABLE X. -continued 







isospin basis 


particle basis 






Expt. 








Coulomb off 


Coulomb on 










-2.51 


-2.48 


-2.27 




-2.42 ± 0.30 


S+p 




4.91 


5.03 


4.56 




3.41 ± 0.30 




at 


0.727 


0.727 


0.834 




0.709 ± 0.001 




n 


-1.20 


-1.29 


-2.82 




-0.783 ± 0.003 






-2.58 


-2.58 


-2.59 




-1.8, -2.0 


Ap 


Ts 


2.83 


2.83 


2.83 




2.8, 5.0 




at 


-1.60 


-1.60 


-1.60 




-1.6, -2.2 




n 


3.00 


3.00 


3.00 




3.3, 3.5 






-2.58 


-2.55 


-2.56 






An 


Ts 


2.83 


2.84 


2.86 








at 


-1.60 


-1.57 


-1.57 








n 


3.00 


3.04 


3.03 










1.41 


1.73 


1.20 








Ts 


-11.0 


-26.6 


-15.0 








at 


1.33 


0.802 


0.914 








n 


-39.0 


-4.90 


-22.7 






TABLE XL 


A summary of the phase-shift behavior 


in each channel related to the A A - 'EN{I = 


1/2) iPi 


- ^Pi channel 


coupling, "disp' 


' in the table denotes a dispersion-like resonance, while "step" denotes a step-like resonance. 








RGM-F 


ESS 


RGM-H 






AN 




step 


step 


disp 




disp 


AN 




disp 


disp 


disp 




disp 


EiV 




<5 < 


5 <0 


<5 ~ ^ 60° 




5 <0 


EJV 3 Pi 




5 <0 


S <0 


6^40° 




5^40° 
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TABLE XII. Contributions to the E p scattering total cross sections from each partial wave at ps = 160 MeV/c, calculated 
in the isospin basis. The results given by fss2 and FSS are listed in the unit of mb. A, B and C denote the scattering processes 
E-p -> E-p, E-p -» E°n and E p — > An, respectively. In the process B, the factor {kf/ki) is not included. 



ABC 





issz 




TOO 


ISSZ 


r bb 




rss2 


r bb 


1 o 1 c 
Do — > OO 


7.3 




6.8 


28.3 


26.5 




3.5 


3.4 


Si — > 6i 


118.2 




117.2 


54.6 


55.8 




60.2 


41.6 


So 3 7-» 

Si — > Di 


0.1 




0.2 


0.3 


0.4 




60.3 


78.0 


3 n I 3 c 


n 1 
U.i 




O.Z 


O.O 


0.4 




U.I 


U.i 


'Po 'Po 


0.5 




0.6 


1.3 


1.4 




0.2 


0.4 


'Pi - 'Pi 


0.1 




0.3 


0.5 


0.8 




2.0 


2.7 


'Pi - ^^Pi 


0.1 




0.3 


0.1 


0.2 




0.5 


2.4 


^Pi ^ 'Pi 


0.1 




0.3 


0.1 


0.2 




8.2 


21.7 


'Pi -» 'Pi 


4.3 




2.4 


4.6 


2.2 




12.2 


15.0 


3 p„ , 3 p 


0.1 




0.1 


0.0 


0.0 




0.1 


U. i 


'P2 'F2 


0.0 




0.0 


0.0 


0.0 




3.3 


4.0 


total (J < 3) 


131.2 




128.5 


90.2 


87.9 




151.1 


169.9 


TABLE XIII. 
Unit of length is 


Scattering length (^4) and effective ranges 
in fm. A'^ and P'^ imply the Coulomb case 


(P) matrices of fss2 for the low-energy E p scattering in 
in the particle basis. 


'So state. 






A 




R 




A" 




R" 


E-p ^ E-p 




0.2644 




2.197 




0.2238 




2.349 


E-p^ E°n 




-1.287 




1.997 




-1.227 




1.845 


E-p An 




-0.1692 




7.231 




-0.1650 




6.683 


E«n ^ E«n 




-0.9617 




2.916 




-1.066 




3.416 


E°n -> An 




0.0887 




-6.165 




0.0755 




-5.016 


An — > An 




-0.0850 




11.39 




-0.0867 




13.50 
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TABLE XIV. Scattering length {A) and effective ranges (7?.) matrices of the reduced expansion for the low-energy E~p 
scattering in the ^Si + state, calculated by fss2 using the particle basis. The ^Di channel of the E~p state is eliminated. 

A'^ and TV^ imply the Coulomb case in the particle basis. Unit of length is in fm^+*'+i in A, A% and in fm^"*-^' in 7^, 11% 
where I and t' are the orbital angular momenta of the initial and final channels, respectively. 





A 


TZ 


^= 


7e= 




-14.91 


2.196 


143.7 


2.527 


T.-p^Si T,°n^Si 


11.25 


-2.594 


-112.7 


-3.274 




-6.500 


9.071 


59.83 


9.204 




5.589 


22.30 


-56.37 


29.97 


T,-p'^Si An^Di 


-2.742 


31.81 


27.36 


44.52 




-7.395 


13.81 


89.46 


14.58 




3.953 


13.92 


-47.82 


15.06 




-4.032 


-153.4 


44.39 


-164.7 




1.979 


-264.2 


-21.55 


-283.8 




-2.166 


28.12 


25.50 


31.76 


TPti^Di An^Si 


1.990 


-88.32 


-23.89 


-102.8 




-0.9492 


-208.6 


11.63 


-236.2 


An^Si ^An^Si 


-1.981 


1664 


22.23 


1821 


An^Si ^An^Di 


0.9603 


2783 


-10.80 


3055 


An^Di -i- An^Di 


-0.4598 


4681 


5.253 


5151 



TABLE XV. S^p zero-energy total cross sections kias (mb- fm~^) for each spin-state (5 = and 1) and ktar for reactions 
B: E~p — > TPn and C: E~p — > An. The divergent factor Co is taten out in the Coulomb case. The inelastic capture ratio vr 
at rest and rp in flight are also given. Two versions of the quark model, ESS and fss2, are used. 







ESS without Coulomb 






fss2 without Coulomb 








B 


C 


rs, Tf 




B 


C 


rs, vf 


tr 


kiCTl 


104 
143 


4.4 
340 


0.959 ] 
0.296 J 


I 0.462 


59.8 
66.7 


5.98 
165 


0.909 ] 
0.289 J 


I 0.444 


ki(TT 


133 


256 


0.342 


65.0 


125 


0.342 






ESS with Coulomb 








fss2 with Coulomb 








B 


C 


rs, Tf 


tr 


B 


C 


rs, rF 


tr 


(Co)~^kiao 


95 
71 


4.4 

168 


0.955 ] 
0.296 J 


I 0.461 


47.7 
38.3 


4.98 
94.6 


0.905 ] 
0.288 J 


I 0.442 


{CoY^kiOT 


77 


127 


0.377 


40.6 


72.2 


0.360 



45 



TABLE XVI. Corrections to the total cross sections in mb; 1) pot: The Coulomb and threshold corrections used for the 
parameter search, 2) fss2: The pion-Coulomb correction (with the correct threshold energies) calculated with AG term included 
in the particle basis. 
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TABLE XVIL A and E s.p. potentials in 
G-matrices in the continuous prescription for 
is also shown for comparison |59|. 


nuclear matter with 
intermediate spectra. 
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TABLE XVIIL 
G-matrices of fss2 


The nuclear-matter density dependence 
in the continuous prescription. The unit 
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TABLE XIX. Decomposition of Sa = -11.1 MeV • fm^ and Ss = -23.3 MeV • fm'' at kp = 1.35 fm"^ into various 
contributions. The model is fss2. Contributions in FSS are also shown in parentheses for comparison. The unit is MeV ■ fm"". 
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TABLE XX. The spatial part of the exchange Born amplitudes defined by Eq. (A. 4) of The polynomial part u{k,q) 
of the two-body force in Eq. (A3) is also shown. The coefficients a, €, A — 2pq/{l — r'^), and the vectors pr = (V/^/^b), 
ar ~ {A/y/2fj,b) are calculated from Eq. (A. 14) of |Q| by setting x = 1 and /i = 3/2 for each interaction type T. The factor 
A is non-zero only for the T — D± types , and e ^ only for the T = S, S' types. The basic spatial functions fr{S) with 
T = C, CD, LS, TD are defined by Eq. (|aTo|). 
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TABLE XXI. The coefficients of^ and the correspondence among il, Q," in the two-bodey force Eq. (Al). The column 
[3 implies the meson types and 7 = {m/2m-ud)- 
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TABLE XXIL The coefficients Cj and Dj in Eq. (06) for the parameterized deuteron wave functions. The model is fss2, 
calculated in the particle basis with the full Coulomb exchange kernel. The number of parameters is n = 11, but the last Cj 
and the last three Dj (the parenthesized values) should be calculated from Eqs. (C.7) and (C.8) of [^. 
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FIG. 1. Threshold relations of the AN - EiV coupled-channel system in the non-relativistic kinematics. 
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FIG. 2. Calculated np phase shifts by fss2 in the isospin basis, compared with the phase-shift analysis SP99 by Arndt et 
al. ^ The dotted curves indicate the result given by FSS. Some empirical inelasticity parameters p of SP99 are also shown for 
Tlab > 300 MeV, in order to give a measure of possible deviations of the phase-shift values in the single-channel calculation. 
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FIG. 3. Decomposition of the ^Pj phase shifts for the np scattering to the central (Sc), LS {5ls) and tensor (St) components. 
The results given by fss2 (solid curves) and FSS (dashed curves) are compared with the decomposition of the empirical phase 
shifts SP99 H. 
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FIG. 4. Calculated np differential cross sections compared with the experiment Calculation is performed using fss2 
with the full pion-Coulomb correction included. 
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FIG. 5. The same as Fig. M but for the np polarization. 
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FIG. 6. The same as Fig.U but for pp differential cross sections. 
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FIG. 7. The same as Fig. ^ but for the pp polarization. 
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FIG. 8. The five invariant amplitudes for the pp scattering at TJab ~ 800 MeV, calculated by fss2 (solid curves), the Paris 
potential (dashed curves) and the empirical phase shifts SP99 |^ (dotted curves). The Coulomb force is included in fss2 
and SP99, but not in the Paris potential. 
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FIG. 9. (a) The momentum dependence of the s.p. potentials Usiqi) predicted by the G-matrix calculation of fss2. 
The QTQ prescription is used for intermediate spectra. The real part Ke!7B(gi) is shown, (b) The same as (a) but for the 
imaginary part 3mf/fl(gi). (c) The same as (a) but in the continuous prescription for intermediate spectra. The nucleon s.p. 
potentials obtained by the p prescription are also shown with respect to the T-matrices of fss2, the Paris potential jsoj and 
the empirical phase shifts SP99 The momentum points selected correspond to Tiab = 100, 200, 400, 800, and 1600 MeV for 
the NN scattering. The partial waves up to J < 8 are included in fss2 and the Paris potential, and J < 7 in SP99. (d) The 
same as (c) but for the imaginary part Ssm)7_B(gi). 
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FIG. 10. The deuteron wave functions predicted by fss2 (solid curves) and by Bonn model-C 
and momentum (lower) representations. 
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FIG. 11. The S- and P-wavc phase shifts of the T.N{I = 3/2) system, predicted by fss2 (sohd curves), (a): '^So, ^Si, (b): 
^Po, ^P2, and (c): ^Pi, ^Pi. Results by FSS (dashed curves) and RGM-H (dotted curves) are also displayed. 
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FIG. 12. (a) The total cross sections and (b) the forward to backward ratio of the differential cross sections (F/B) for the Ap 
elastic scattering, predicted b y fs s2 in the isospin basis. The results by FSS are also shown in dashed curves. The experimental 
data are cited from Refs. 




FIG. 13. The AN ^Si and ^5*0 phase shifts calculated by fss2 (left) and by FSS (right) in the isospin basis. 
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FIG. 14. The S'-matrix of the AA'' - T,N{I = 1/2) ^Pi - ^Pi coupled-channel system, predicted by our three models fss2 
(upper), FSS (middle) and RGM-H (lower) in the isospin-basis. The diagonal phase shifts Su (left-hand side), and the reflection 
and transmission coofEcicnts rjn (right-hand side) defined by Sij = rjije^^^'i are displayed. The channels are specified by 1: AN 
iPi, 2: AN ^Pi, 3: T.N and 4: T.N ^Pi. 
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FIG. 15. (a) The AN elastic total cross sections and (b) AN — » UN reaction total cross sections in the threshold region, 
predicted by fss2 with the full pion-Coulomb correction. In (a), the AN elastic total cross section in the isospin basis (without 
the Coulomb force) is also shown by the dotted curve. 
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FIG. 16. The same as Fig.|i^ but for the ^^o phase shift and for the ^-matrix of the KN - T.N{I = 1/2) ^Si - ^Di 
coupled-channel state. In the coupled-channel system, the channels are specified by 1: AA*' ^S\, 2: AA^ '^Di, 3: EA'' Si and 4: 
T.N ^Di. The upper figures display the result given by fss2, while the lower by FSS. 
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FIG. 17. Argand diagram of the S'-matrix element S33 for the AA*' - EA(J = 1/2) - ^Di coupled-channel system. The 
channel i = 3 corresponds to the 'SN{I = 1/2) "^Si channel. The range of the incident momentum is < pe < 1000 MeV/c. 
Results given by fss2 (solid curve), FSS (dashed curve) and RGM-H (dotted curve) are displayed. 
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FIG. 18. (a) The low-energy E~p, YPn and An phase shifts for the ^So channel, predicted by fss2, as a function of the 
incident momentum ps. The solid curves include the full pion-Coulomb correction i n the particle basis, while the dashed curves 
are obtained in the effective range approximation using the parameters in Table Kill, (b) E~p reflection and transmission 
coefficients r]i\ for the ^So channel, predicted by fss2. The final channel i is specified by 1: E~p, 2: YPn and 3: An. The solid 
curves denote the full calculation, while the dashed curves the effective range approximation. 
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FIG. 19. (a) E p, E°n and An phase shifts for the S\ + '^Di state, predicted by fss2, as a function of the incident 



momentum ps. E" 
in the solid curves 
(b)E 
by 1: E-p Sx, 2 



p D\ phase shift is not shown, since it is smaller than 0.1°. The full pion-Coulomb correction is inc luded 

The dashed curves are obtained in the effective range approximation using the parameters in Table 
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p reflection and transmission coefficients r]i\ for the ^5*1 -I- '^-Di channel, predicted by fss2. The final channel i is specified 



E-p •'Di, 3: E" 



?i, 4: E n -Di, 5: An Si and 6: An D\. 7721 and 7^41 are not shown, since they are 



very small. The solid curves include the fuU pion-Coulomb correction, while the dashed curves are obtained in the effective 
range approximation. 
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FIG. 20. E p inelastic capture ratio in flight rp, predicted by fss2, as a function of tfie incident momentum ps. Calculation 
is made in the particle basis with (solid curve) and without (dashed curve) the Coulomb force. 
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FIG. 21. Calculated low-energy of E+p and E~p scattering total cross sections by fss2 (solid curves), compared with 
experimental data: (a) Tj'^p elastic, (b) E~p elastic, (c) E~p — > TPn charge-exchange, (d) E~p An reaction cross sections. 
Predictions in the particle basis without the Coulomb force (dashed curves), and those in the isospin basis (dotted curves) are 
also shown. The experimental data are taken from [p3| for (a) and (b), and from fc2| for (c) and (d). 
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FIG. 22. Calculated E p and S~p differential cross sections by fss2, compared with the experimental angular distributions: 
(a) E"'"p elastic scattering at ps = 170 MeV/c, (b) the same as (a) but at ps ~ 450 MeV/c, (c) S~p elastic scattering at 
Pe = 160 MeV/c, (d) The same as (c) but at ps = 550 MeV/c, (e) E~p ^ An differential cross sections at ps = 160 MeV/c. 
The experimental data are taken from [B2| for (a), [M for (b), for (c) and (e), and |B5| for (d). 
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FIG. 23. Calculated A'^A'^ and YN total cross sections by fss2, compared with the experimental data. The solid curves 
denote the full calculation and the dotted curves the calculation in the isospin basis. The experimental data are taken from ^] 
for NN, H, @, H for Ap, H for E+p, and H, g], H for E'p scattering. 



67 




fss2 (cont.) 



1.2 



1.4 



1.6 1. 
kp [fm-^] 



FIG. 24. Nuclear matter saturation curves obtained for fss2 and FSS, together with the results of the Paris potential ]5€| ] 
and the Bonn model-B (Bonn-B) potential js^. The choice of the intermediate spectra is specified by "QTQ" and "cont.", 
respectively. The result for the Bonn-B potential in the continuous choice is taken from the non-relativistic calculation in . 
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FIG. 25. (a) The nucleon s.p. potential I7jv(gi) in nuclear matter in the continuous choice for intermediate spectra. 
Predictions by fss2 for three densities p = 0.5 po, 0.7 po and po are shown. Here the normal density po corresponds to 
kp = 1.35 fm-^ The dashed curve is the result achieved by Schulze et al. |^ with the Nijmegen soft-core A''A'^ potential 
NSC89 jssj. (b) The same as (a) but for the A s.p. potential Uk{(1i)- (c) The same as (a) but for the E s.p. potential UT.{q_i)- 
(d) The same as (a) but for the imaginary part of the A and E s.p. potentials I7fl(gi) predicted by fss2. 
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